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SECTION  I 
INTRODUCTION 

Propagation  of  infrared  laser  radiation  through  the  atmosphere  is  to 
a greater  or  lesser  extent  hindered  by  the  presence  of  water  vapor.  Spe- 
cifically, molecular  absorption  on  the  various  rotational- vibrational 
levels  of  water  give  rise  to  significant  atmospheric  attenuation  of  high- 
intensity  laser  radiation  at  CO2,  CO,  and  HF  frequencies  and  to  a lesser 
extent  at  DF  frequencies.  A detailed  examination  of  the  molecular  struc- 
ture of  water  vapor  indicates  that  the  following  vibrational  bands  arc  of 
importance:  the  ground  state,  (010),  1.100),  (001),  (020)  as  well  as 
mixed  states  (110)  and  lOll).  Different  regions  of  the  water 
vapor  absorption  spectrum  involve  different  vibrational  bands. 

Furthermore,  due  to  the  extreme  asymmetry  of  this  molecule  as  well  as  its 
large  rotation  constants,  the  level  structure  is  sufficiently  complex  that 
often  transitions  from  the  same  band  cover  different  regions  of  the  spec- 
trum. For  example,  ground  to  transitions  arise  in  both  the  CO2  and 
CO  regions  of  tiie  spectrum,  and  this  in  conjunction  with  the  extreme 
centrifugal  distortion  that  occurs  in  water  demands  tliat  cne  treat  tlie 
molecular  rotations  with  great  care.  Furthermore,  as  noted  abo\e,  the 
different  regions  of  the  spectrum  involve  different  vibrational  bands, 
and  the  different  vibrational  states  interact.  (See  Tables  1 and  2.) 
Consequently,  a multimode  model  of  the  water  molecule  is  required 
for  a proper  description  of  those  regions  of  the  water  vapor  spectru.m 
tnat  are  of  interest.  In  this  report  we  describe  three  different  ap- 
proaches used  in  developing  a multimode  model  of  water  vapor  along  witi'. 
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Table  1 


Spectrum 

of  Water  Vapor  from  2 

to  10.6  am 

Frequency 

Wavelength 

(um) 

Description 

C0-, 

10.6 

kinetic  cooling,  V-V  rates 

of  N'.,  off  H.,0  water  vapor 

continuum 

CO 

4. 8-5.0 

water  vapor  continuum 

DF 

5-4 

water  vapor  continuum  and 

HDO 

HF 

2-5 

water  vapor  lines 

Table  2 

Vibrational  Characteristics  of  CO  and  CO? 


CO 


CO2 


Numerous  hot  band  transitions 
at  5 jjm:  v?  -►  Zv? 

Many  cross-band  transitions 


involving 

the  hot  band: 

V2 

- V3 

-*  'J3 

Ground 

to 

vj  transitions 

i 

. V2 

-*  Vl 

1 '•'2 

-*  V3 

At  4.5 

1 

\0 

— v; 

1 

1 

'■*2 

“ -'->2 

Hot  band  transitions  at  9.28  am: 
Ground  to  v?  dominates  10.6  am 

Ground  to  v?  dominates  9.23  am 

Pure  rotational  transitions  at 
9.2s  -m: 

11-5  - 1^3 

15-’i  - 14-3 

i-i-io  - 


15-3 


relevant  results.  The  first  involves  the  use  of  a single  coordinate  tha 
is  so  designed  to  simulate  all  three  vibrational  degrees  of  freedom  and 
is  applied  to  the  ground  and  (010)  vibrational  states.  'These  states 
play  a significant  role  for  CO2.  CO,  and  DF  frequencies.  The  second 
approach  involves  treating  all  three  vibrational  coordinates  as  well  as 
their  interactions.  Finally,  a third  approach  that  is  motivated  by  a 
vibrational  analysis  is  presented.  This  technique  is  phenomenological 
in  nature  and  offers  e.xcellent  hope  for  a global  (1  to  20  um)  model  of 
water  vapor  that  is  highly  accurate  and  reliable. 

In  Section  II  we  present  our  results  for  the  correlated  bender 
approach  to  water  and  discuss  the  physical  content  of  our  model.  In 
Section  III  we  discuss  our  predictive  phenomenology  for  a multimode 
treatment  of  water  and  display  preliminary  results.  Finally,  in  Section 
ri  we  present  a multimode  treatment  of  the  water  molecule,  describe  the 
computer  programs,  outline  the  physical  content  of  the  model  and  exhibit 


its  results. 


SECTION  II 

CORRELATED  BENDER  MODEL  OF  WATER  VAPOR 


( 
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In  this  section  we  present  a model  for  treating  the  rotational- 
vibrational  structure  and  spectra  tliat  arise  from  bond  angle  fluctuations 
in  molecular  systems  and  appl\'  it  to  the  specific  case  of  water.  Our 
model  is  based  on  the  motion  of  the  bond  angle  and  the  correlation  of  this 
motion  with  tiie  stretch  modes  and  involves  the  description  of  a single 
intsrnuclear  coordinate  (but  is  multimode  ' nature''. 

Specifically,  we  present  a model  of  the  ater  molecule  that  is  based 
on  tne  relative  motion  of  the  two  protons,  with  respect  to  each  other,  and 
tiie  coupling  of  tiiis  motion  to  the  molecular  rotation.  In  this  model,  the 
classical  trajector>'  of  the  two  protons  is  a straight  path.  T'nis  implies 
that  the  OH  bond  length  is  no  longer  fi.xed  at  some  average  value,  but  is 
allowed  to  fluctuate  along  v^ith  the  bond  angle.  Since  we  are  describing 
tnis  motion  by  means  of  a single  coordinate,  it  is  clear  that  such  a model 
contains  correlations  between  the  OH  bond  length  and  the  bond  angle  fluc- 
tuations. The  most  interesting  feature  of  this  model,  which  we  shall  refer 
to  as  a correlated  bender  model,  is  tiiat  it  ;.'ields  a more  accurate  descrip- 
tion of  the  rotational- vibrational  level  structure  of  the  ground  and  (010) 
vibrational  states  in  the  seni-c  that  its  predicted  energies  and  molecular 
geometry  are  in  better  agreement  with  experiment  than  previous  work.  Typ- 
ical values  are  within  O.Sii  for  all  rotational  states  th.rough  J = 10  and 
tiirougii  J = 15.  In  constructing  this  model,  the  only  e.xperimental  data 
used  were  the  first  three  (.J  = 0)  bending  mode  vibrational  energy  level 
spacings  as  well  as  the  A and  3 rotation  constants  associated  with  the 


4 
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ground  vibrational  state.  Note  that  the  C rotation  constant  of  the  ground 
state  as  well  as  all  of  the  rotation  constants  of  the  (010)  vibrational 
state  are  predicted  quantities.  It  is  found  that  the  correlated  bender 
predicts  a ground  state  C-rotation  constant  that  lies  within  0.01  cm*^  of 
the  measured  value.  Furthermore,  the  e.xcited  state  rotation  constants  are 
found  to  lie  within  l°o  of  the  experimental  values. 

In  Section  IIA  we  present  our  model  of  water  vapor  and  discuss  its 
physical  characteristics.  Specifically,  we  begin  with  the  Darling-Dennison 
rotational -vibrational  Hamiltonian  for  a polyatomic  molecule  (ref.  1)  and 
in  a series  of  steps  reduce  it  to  two  different  forms  that  are  appropriate 
for  treating  our  two  models  of  water.  Following  this  we  present  a detailed 
description  of  the  vibrational  structure  of  H2O  that  our  model  predicts. 
Specifically,  we  present  tiie  potential  energy  surface,  the  characteristic 
structure  of  the  proton-proton  motion,  and  the  rotation-constant  matrices 
associated  with  this  model. 

In  Section  IIB  we  consider  the  rotational  level  structure  within  the 
ground  and  (010)  vibrational  bands.  Specifically,  we  obtain  the  vibra- 
tional-rotational energies  and  wave  functions  using  a modified  theory  of 
the  asymmetric  rotor.  In  particular,  we  diagonalite  the  vibrational- 
rotational  energies  toget’p.er,  not  individual!;/  as  is  usually  done.  Tliis 
combined  diagor.alitation  is  particularly  important  for  K2O  inasmuch  as 
the  rotational- vibrational  interaction  is  extremely  strong  and  consequently 
plays  an  important  role  in  tiie  rotational -vibrational  level  structure. 

We  present  tables  of  predicted  versus  experimental  values  of  the  various 
rotational  states  of  ilgO.  To  determine  how  much  of  the  rotational- 
vibrational  interaction  is  accounted  for  in  our  models,  we  also  include 


r 


in  these  tables  the  values  of  the  rigid  rotor  energies.  These  were 
obtained  by  replacing  all  off-diagonal  elements  of  the  rotation  matrices 
with  :ero  (ref.  2).  .A  comparison  of  the  various  energies  demonstrates 
that  for  the  lOi^  state,  about  90°i  of  the  rotational-vibrational  inter- 
action is  accounted  for  in  the  correlated  bender  model. 

In  Section  IIC  we  discuss  the  symmetric  isotopes  D2O  and  TiO.  This 
enables  us  to  e.xamine  certain  features  of  the  correlated  bender  model  that 
are  mass  dependent,  specifically,  the  degree  of  centrifugal  distortion, 
the  molecular  structure,  and  the  correlation  between  the  bond  angle  fluc- 
tuations and  the  OH  stretch  motion. 


.■\ . Foundations  for  Quantim  .-VnalNsis 

In  this  section  wo  present  our  model  of  t!ie  water  molecule  that  is 
based  on  the  fluctuations  of  a single  coordinate.  For  convenier.ee , we 
have  subdivided  this  section  into  two  parts.  In  Section  II.-\-l  we  con- 
struct tlie  quantum  mechanical  Hamiltonian  that  we  use  to  model  the 
rotational-vibrational  structure  and  spectrum  of  the  water  molecule.  In 
Section  IIB-1  we  use  the  Hamiltonian  to  derive  various  properties  of  H^O 
that  depend  on  tire  intranuclear  motion. 

1.  Rotational-Vibrational  Hamiltonian 


Our  starting  point  is  the  Darling-L'ennison  (ref.  1)  rotational- 
, ibrational  Hamiltonian  for  a nonlinear  pol>Mtomic  molecule 
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a y.  ao 


k ' p,  u''"  P.  a''-*  * V 
k ^ 


CD 


0 


In  liquation  (.1),  \ is  the  molecular  potential  energy  surface,  ^hich  is  a func- 
tion of  all  the  normal  coordinates  = 1.  Pj,  is  the  correspon- 

ding conjugate  momenta,  I!  is  the  ath  component  of  tiie  total  rotational- 
angular  momentum  of  the  molecule,  and  is  the  ath  component  of  the 
vibrational-angular  momentum,  is  tlie  [a,Sj  component  of  the  effec- 

tive reciprocal  moment  of  inertia  tensor,  and  u is  its  determinant.  These 
quantities  are  determined  by  the  molecular  configuration  and  therefore 
are  functions  of  tiie  normal  coordinates. 

Watson  (ref.  5;  has  demonstrated  that  for  nonlinear  molecules  the 
Darling-Dennison  Hamiltonian  can  be  cast  in  the  following  very  useful  form 


H = 4 I c 


at 


U “ ' 

a 


i V p 2 

- k 


C) 


where  U is  given  b\- 


U 


(3) 


It  follows  from  liquation  i.i)  that  U depends  onl>-  on  the  nuclear  masses,  the 
molecular  geometry,  and  tiie  normal  coordinates.  Hence,  it  is  independent 
of  tiie  various  momenta  tiiat  appear  in  liquation  ill,  and  as  noted  by  Watson, 
can  be  regarded  as  an  effective  mass-dependent  potential  tliat  we  formally 
include  in  V.  We  note  tiiat,  when  a molecule  assumes  a linear  configura- 
tion, the  nuclear  moment  of  inertia  associated  with  that  direction  vanislies 
and  as  a consequence. 


LI  - 


(4) 


The  presence  of  an  infinitely  attractive  term  in  a molecular  Hamiltonian 


is  disturbing  and  on  physical  grounds  should  not  be  present.  The  resolution 


of  this  difficulty  can  be  obtained  by  noting  that  for  linear  molecules 
tile  rotational-vibrational  Hamiltonian  can  be  written  as  (ref.  4) 


H = 


X 


\x  " i "■^^x)(^x-\x)-  ^ I 

(^v-y  ^ 3 ^ H ^ 


(5) 


where  the  quantities  and  are  defined  in  Reference  4 and  do  not  directly 
concern  us  here.  The  important  feature  of  Equation  (5)  is  the  absence  of  the 
effective  potential  U.  We  note  that  Equations  (21  and  (5)  differ  considerably 
in  form,  particularly  with  regard  to  the  rotational  structure  [n_  is  not 
present  in  (5)]  and  the  fact  that  linear  molecules  have  an  additional 
vibrational  degree  with  respect  to  nonlinear  molecules.  These  features 
are  well  known;  our  only  purpose  in  discussing  them  is  to  point  out  that 
tne  molecular  Hamiltonian  is  well  behaved  in  the  sense  tliat  it  does  not 
contain  terms  that  diverge  to  minus  infinity  (ref.  4) . In  the  remainder 
of  this  paper  we  siiali  group  the  effective  potential  U with  the  actual 
potential  V and  refer  to  the  sum,  U + V,  as  V.  It  is  this  quantity  that 
will  be  obtained  by  a least-squares  fitting  technique.  We  note,  as  a 
consequence,  that  tlio  potential  energy  that  we  eventually  obtain  is  mass 
dependent.  We  can  investigate  this  mass  dependence  by  examining  the  iso- 
topic species  D2O  and  T2O.  This  is  done  in  Section  IIC.  Thus,  Equation 
[2'.  now  reads 


t.iat 


H 


(n  -TT  )u  .(2.--J 


For  a symmetrical  molecule  such 
tne  reciprocal  moment  of  inertia 


as  H2O,  it  is  easy 
tensor  is  diagonal 


Ca) 


to  demonstrate 
thus  the 


S 


molecular  rotational- vibrational  Hajniltonian  reduces  to 


H = y y (n  -TT  )u  (n  -TT ) 4 y P,  ^ 

2 ^ a a aa  a a ^ k 


+ V 


(6) 


Furthermore,  we  shall  neglect  both  Coriolis  forces  as  well  as  vibrational- 
angular  momentum.  This  appro.ximation  is  justified  on  the  grounds  that  the 
bulk  of  the  rotational-vibrational  interaction  in  the  (000)  and  (010) 
bands  arises  fi'om  centrifugal  distortion.  Thus,  our  molecular  Hamiltonian 
reduces  to 

1 r-  ^ 1 „ „ 

(7) 


H = 4 I (7  Vu  ) * k I PJ  * V 


k 

.\e.\t,  we  consider  the  kinetic  energy  associated  with  the  bending 
mode  vibration,  i.e.,  fluctuations  in  the  bond  angle  0.  We  begin  by  noting 
that  tile  moment  of  inertia  1(6)  associated  witii  the  bending  mode  is 


1(6) 


1 7 r 


ito  m|,(l  + cose) 


/M. 


(8) 


where  r is  the  OH  bond  length,  mL,,  m„,  and  M,,  tiie  atomic  hvdrogen,  ox''-- 

n U ^2^ 

gen,  and  molecular  water  masses,  respective!)’ . We  note  tiiat  1(e)  depends 
only  weakly  on  the  bond  angle  e and  in  Section  I IB,  we  demonstrate  that 
tiie  contribution  to  the  vibrational  energy  that  arises  from  fluctuations 
in  I(S)  amounts  to  less  than  O.OSo.  .Now,  the  Born-Oppenheimer  theorem 
(ref.  5)  states  that  the  coupling  between  the  electron  and  nuclear  motion 
is  on  tile  order  of  (m^/Mj’^,  where  M is  some  suitable  nuclear  mass.  Tiius , 
errors  encountered  in  neglecting  fluctuations  of  1(B)  appear  to  be  of  the 
same  order  of  magnitude  as  corrections  that  arise  from  deviations  from  an 
adiabatic  model.  Furtiiermore , as  we  shall  demonstrate  in  Section  IIB,  this 
energy  is  an  order  of  magnitude  smaller  tlian  any  other  that  arises  from 
approximations  made  in  this  paper.  Consequently,  we  shall  ignore  the 


angular  dependence  of  1(9)  and  use  for  the  kinetic  energy  of  the  bending 
mode 


T 


h2 

21  (9o) 


3^ 

'W 


where  lo  is  the  moment  of  inertia  evaluated  at  a fixed  angle  9o  (see 
below) . 

Next,  we  note  that  for  situations  involving  large  amplitude  changes 
in  the  bond  angle,  as  occurs  in  water,  it  is  more  convenient  to  work  with 
a length  than  an  angle.  Thus,  we  define 

X = 2r  sin9/2,  (9) 

which  varies  from  cero  to  2r  as  6 varies  from  zero  to  m . Defining  the 
dimensionless  coordinate  q as  x/r,  we  find  that  the  kinetic  energy  asso- 
ciated ivdth  fluctuations  of  the  bond  angle  is 


Fixing  the  OH  bond  length  at  r,  we  obtain  the  following  Hamiltonian 

2 


a aa  ^ ^ 


where  V],  (q)  will  be  taken  to  be  a power  series  of  the  form 

''■l(q)  = I a^[q'qo]'"‘ 

n=2 


(12) 


Note  that  Equation  (12)  corresponds  to  an  expansion  in  Legendre  polynomials, 
P^(cos9),  which  form  a complete  set  where  q;  = 2sin9g/2. 

B.  Vibrational  Structure 

In  our  correlated  bender  model,  we  wish  to  constrain  the  classical 
trajectory  in  such  a way  that  the  two  protons  will  travel  in  straight 


10 


paths.  Now,  the  variation  of  the  kinetic  energy  operator  with  q as  well 
as  the  presence  of  the  first  derivative  term  q(d/dq)  in  Equation  (10)  reflect 


'hat  the  classical  trajectory  of  the  two  protons  is  along  the 
cle  of  radius  r.  Thus,  if  we  neglect  the  first  derivative 
ter:ii  .d  i.r  ' q to  a value  qo  (discussed  below),  then  classically  we  are 
describing  the  motion  of  two  particles  of  mass  u 


”’o 

+ m,. 
n 

qo/2) 

2' 

1 _ 

1 - 1 

[qo/2)-- 

1 . 598  m 


(15) 


traveling  along  a straight  path.  The  molecular  Hamiltonian  H for  this 
model  is 


H 


ct 

ct  cca 


1(2  d 
ur-  dq^ 


V(q) 


(14) 


where 


V(q) 


S ^n 


(15) 


qQ  = 2 sin9Q/2,  r and  a^  are  so  chosen  that  the  eigenvalues  and  eigen- 


vectors of  H will  least-squares  fit  tiie  A and  B rotation  constants  of 
the  ground  vibrational  state  as  well  as  t.he  spacings  between  ground  and 
first  three  e.Kcited  bending  mode  states.  This  completel\-  defines  our 
model . 


At  first  glance,  one  might  suppose  that  our  model  as  mani- 
fested in  Equation  (14)  represents  an  appro.xi.mati  on . Till s , however , is  not 
the  case.  Specifically,  we  are  now  including  correlations  between  fluc- 
tuations of  the  OH  bond  length  and  the  bond  angle.  This  follows  immedi- 
ately if  one  considers  the  structure  of  the  motion  of  the  protons.  Spe- 
cifically, for  configurations  in  which  the  protons  extend  beyond  the 


11 


equilibrium  distance  qg  the  OH  bond  distance  increases  whereas  if  q < qg, 
this  distance  decreases.  Thus,  the  fluctuations  in  the  bond  angle  and 
bond  length  are  strongly  correlated  in  our  model.  Note,  however,  that  we 
do  not  take  into  account  the  complete  motion  of  the  OH  bond  length.  Thus, 
this  correlated  bender  model  should  be  regarded  as  an  alternative  descrip- 
tion of  the  water  molecule.  This  interpretation  is  clearly  justified  by 
the  fact  that  our  model  renders  a highly  accurate  description  of  the  H2O 
molecule . 


1.  Vibration  Structure  of  the  Water  Molecule 

In  this  section  we  present  the  vibrational  properties  of  the  water 
molecule  that  are  predicted  by  our  model.  Specifically,  we  displag'  tables 
of  (1)  the  potential  energy  constants,  (2)  e.xpectation  values  of  various 
moments  of  the  vibrational  coordinate,  (3)  the  .A.,  B,  and  C rotation  con- 
stant matrices,  and  (4)  the  geometric  structure  of  the  molecule.  For 
convenience,  these  properties  are  displayed  in  separate  sections.  Before 
doing  so,  we  first  describe  our  least-squares  (self-consistent)  fitting 
technique.  Our  basic  approach  is  to  first  determine  the  J = 0 bending 
mode  vibration  eigenfunctions  and  eigenvalues.  This  is  done  by  iterating 
the  following  steps  self-consistently. 

The  vibrational  eigenvalues  (E^)  and  eigenfunctions  deter- 

mined by  first  solving  the  SchrSdinger  equation 


)-  21g 


_d^ 

dq2 


u^iq) 


" En^n(q) 


(lb) 
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Taking  these  eigenfunctions,  we  then  evaluate  the  A and  B rotation 
constants  in  the  ground  and  (0,1,0)  states.  ',Ve  then  determine  the  opti- 
mum pQ,  mean  OH  bond  length,  and  a^  by  demanding  that  the  eigen- 
values and  eigenfunctions  of  (16)  and  (17)  least-squares  fit  the  above- 
mentioned  energy  spacings  and  rotational  constants  (refs.  1,2,6). 

Potential  Energy  Function.  In  Table  3 we  display  the  various  poten- 
tial energy  constants  for  our  model.  N'ote  that  the  bending  mode  (de- 
scribed by  0)  does  not  undergo  simple  harmonic  motion  since  the  bond  angle 
3 is  related  to  q by  means  of 

0=2  sin*^  (q/-) 


Table  5 

Potential  Energy  Constants  Display  [cm'^j 


Para.meter 

Correlated  Bender 

45510 

aa 

0 

ac 

-30833 

3-5 

0 

-25665 

a? 

0 

aa 

1S20 

We  note  that  including  terms  of  the  form  a9(q-qo)^  and  a; 3 (q-qo ) * '^  did 
not  significantly  improve  the  fit.  .At  first  glance,  one  might  suppose 
that  a power  series  e.xpansion  in  (q-qo)  would  be  sufficient  to  fit  three 
level  spacings.  However,  our  scheme  requires  that  the  ground  state  A and 


15 


r 


B rotation  constants  be  accounted  for  as  well,  and  this  affects  the  fitting 
to  the  vibrational  energy  level  structure,  particularly  since  the  greatest 
deviation  appears  in  the  ground  to  first  transition  frequency. 

Proton-Proton  Motion.  .Next,  we  examine  the  ground  state  expectation 
values  of  the  even  moments  of  (q-qo)  (Table  4) . The  first  of  these  quan- 
tities is  related  to  fluctuations  in  the  dipole  moment  and  all  play  a 
role  in  the  rotation  constants. 


Table  4 

Ground  State  Even  .Moments  of  (q-qo) 


,^Ioment 

Correlated 

Bender 

(q-qo)" 

0.910  X 

10"^ 

(q-qo)"* 

0.250  X 

10-2 

(q-qo)° 

0.114  X 

10-^* 

(q-qo)^ 

0.721  X 

10-^ 

(q-qo)"^ 

0.578  X 

10"" 

Rotational  Constants.  It  is  well  known  that  the  rotational  level 
structure  within  a given  band  is  determined  by  the  rotational  constant 
matrix.  In  Table  5 we  display  the  A,  B,  and  C rotational  constant  matrices 
that  are  obtained  from  a correlated  bender  model  of  H^O.  Finally,  in 
Table  6 we  display  the  measured  experimental  values  for  the  diagonal  ele- 
ments of  the  A,  B and  C rotational  constants  (ref.  6) . 

A very  interesting  feature  of  this  state  is  that  both  the  A and  B 
rotational  constants  increase  relative  to  the  ground  state.  Now,  on  naive 
grounds  one  might  e.xpect  that  as  the  molecule  opens  up  (i.e.,  the  bond 
angle  increases)  the  B-rotational  constant  should  decrease.  The  fact 
that  it  does  not  is  indicative  of  a correlation  between  the  ground  state 
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Table  5 


Rotational  Constant  Matrix  for  Correlated  Bender 


A 

Rotation 

Constants  = 

0 

27.876 

-6 . 558 

2.504 

1 

-6 . 358 

31.470 

-11.167 

2.504 

-11.167 

56 . 394 

B 

Rotation 

Constants  = 

0 

14.510 

1.781 

0.227 

1 

1.7S1 

14.837 

2.632 

o 

0.227 

2.652 

15 . 156 

C 

Rotation 

Constants  = 

0 

9.29:^ 

0.084 

-0.009 

1 

0.084 

9.284 

0. 145 

-0.009 

0.143 

9.214 

Table  o 

Experimental  Rotation  Constants 

A Rotation  Constant  B Rotation  Constant  C Rotation  Constant 

:7.S7o  --  14. 5i:  --  9.285 

31.12  --  14.66  --  9.15 
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OH  stretch  mode  and  the  first  excited  bending  mode.  This  interpretation 
is  underscored  by  our  calculation.  Specifically,  the  correlated  bender, 
which  does  include  this  type  of  correlation,  predicts  a larger  B rota- 
tional constant,  which  is  in  agreement  with  experiment. 

Next,  we  note  that  the  strength  of  the  rotational -vibrational 
interaction  is  set  by  the  siie  of  the  off-diagonal  matrix  elements  of  the 
rotational  constants.  For  example,  from  the  correlated  bender  model  we 
have  for  the  .A- rotational  constant 


and 


! i> 

= 0.22 

<1 

'a!  2> 

W 

:0> 

<1 

1a1  1> 

the 

B 

co: 

ns  tant , 

B 

= 0.12 

<1 

|b|2> 

B 

|o> 

<1 

iBll> 

<0 

C 

il" 

= 0.01 

liJ 

C|2> 

<0 

C 

i0> 

<1 

C i 1> 

= 0.56 


= 0.  IS 


= 0.02 


for  the  C constant.  .Note  the  dramatic  increase  in  these  ratios  as  we  go 
from  the  ground  to  the  first  excited  state,  which  is  to  be  expected  as 
the  molecule  opens  up  and  becomes  less  rigid. 

.Molecular  Geometry.  Finally,  we  exa.mine  the  geometric  structure 
that  our  models  of  the  water  molecule  imply.  We  recall  that  the  struc- 
ture parameters  are  so  chosen  that  the  model  reproduces  the  experimentally 
obtained  ground  state  A and  B rotational  constants,  '.'.'e  take  t!ie  static 
bond  angle  (bond  angle  in  the  absence  of  cero  point  vibrations)  and 
ground  state  bond  angle  from  the  relations 


lb 


2 sin 


(qo/^) 


iQ  = 


(i:) 


<0|9|0>  = 2<0|sin-^  (13) 

In  Table  ; we  present  predicted  versus  measured  values  for  the  bond 
length,  bond  angle,  and  63  for  the  model , An  examination  of  the  data  re- 
veals that  the  correlated  bender  yields  a static  bond  angle  that  lies 
within  O.OSijof  tno  experimental  value.  Furthermore,  the  equilibrium  bond 
angle  is  within  0.09'’>of  the  measured  value. 


Table  7 


Geometric 

Structure  of  Water 

Parameter 

Observed 

Correlated  Bender 

Deviation 

(‘’■-) 

r 

0.972  A 

0.960  X 

0.51 

<7> 

104.50° 

104.41° 

0.09 

eo 

105.9° 

105.95° 

0.05 

We  can  also  demonstrate  that  fluctuations  of  the  moment  of  inertia, 
it-),  associated  with  the  bending  motion  are  rot  significant.  T!ie  site 
of  this  effect  is  set  by  the  ratio 


<Ii.'3)>  - 1(9:}  j 


S ^ 10--^  . 


This  corresponds  to  an  energy  of  0.52  cm'l  (ti-.e  ground  state  energ.v  of 
H;0  is  about  850  cm'^)  and  does  not  represent  a sigi.ificant  effect  in  the 
rotational-vibrational  structure  of  the  molecule.  Final  I;.’  the  hydrogen- 
hydrogen  separation  for  the  correlated  bender  is  1.55  X. 
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C.  Rotational-Vibrational  Structure  of  Water  | 

As  we  discussed  in  Sections  IIA  and  IIB,  the  rotational  structure  of  ! 

! 

the  water  molecule  is  dominated  by  centrifugal  distortion.  Since  our 
model  directly  incorporates  this  feature  of  HoO,  an  e.xcellent  test  of  its 
accuracy  is  the  rotational-vibrational  energy  levels  themselves.  In  this 
section,  we  present  tlieoretical  calculations  through  J = 10  of  tiie  ground 
and  first  e.xcited  bending  mode  states  using  our  model  of  water.  In  gen- 
eral, we  find  that  the  correlated  bonder  displays  a reliable  picture  of 
water  as  its  predicted  energy  levels  lie  quite  close  to  e.xperimental  find- 
ings. Furthermore,  to  measure  the  significance  of  the  off-diagonal  rota- 
tional constant  matrix  elements,  we  also  present  predicted  rotational 
energies  for  a model  in  which  these  matrix  elements  have  been  set  equal 
to  tero.  We  note  that  these  matrix  elements  represent  the  coupling  i 

between  the  various  vibrational  states  and  is  a measure  of  the  imuortance  | 

of  centrifugal  distortion  on  the  structure  of  the  molecule.  These 
energies  can  be  regarded  as  rigid  rotor  energies  despite  the  fact  that 
the  rotational  constants  are  given  by  the  reciprocal  moment  of  inertia 
tensor  averaged  over  the  ground  or  first  vibrational  state. 

In  rable  S we  display  the  predicted  rotational-vibrational  level 
structure  in  the  ground  vibrational  state  through  J = 10  (ref.  “) . As 
we  are  specifically  interested  in  the  effects  of  centrifugal  distortion 
on  the  structure  of  the  molecule,  we  present  only  those  states  with  t > 0. 

It  is  wei’  known  that  such  states  involve  rotatic.;  about  the  small  mo- 
ments of  inertia  [A  and  B) . On  physical  grounds  it  is  clear  that  the 
coupling  between  molecular  rotation  and  vibration  is  strongest  for  these 
states.  This  feature  is  reflected  in  the  predicted  lei'el  structure. 
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Co:nparison 

of  Rotati 

onal-Vibrat ional  Energy  Levels  in 

COOO,^  St 

J 

Rigid 

Correlated 

T 

Observed 

Rocor 

5 Dsv . 

bender 

V Dev. 

ii 

42.37 

42.58 

o.co 

42.37 

0.00 

2o 

95.13 

95.20 

0.02 

95.21 

0.05 

22 

136.16 

156.55 

0.28 

136.18 

0.01 

206.30 

206.69 

0 . 33 

206.48 

0.09 

^3 

285.41 

287.26 

0.65 

285.52 

0.04 

•’o 

515.78 

516.45 

0.12 

316.14 

0.11 

42 

333.84 

585.79 

0.51 

584.56 

0.14 

44 

488.13 

494.49 

1.30 

488.43 

0.07 

So 

503.97 

506 . 17 

0.44 

505 .02 

0.21 

= 3 

610.34 

616.49 

l.CO 

611.58 

0.20 

5 = 

742.07 

757.39 

2.06 

745.04 

0.13 

Oo 

661 .55 

664 . 80 

0.49 

665.14 

0.24 

62 

757 . 7S 

764.17 

0.84 

760.03 

0.50 

64 

8SS.63 

905.43 

1 . O’ 

891.24 

0.29 

66 

1045.06 

1076.05 

2.96 

1047.30 

0.21 

'0 

927.75 

935.07 

0.79 

951.50 

0.33 

23 

1059.85 

1074.59 

1.59 

1064.10 

0.40 

’5 

1210.19 

1240.28 

1 « - 

1221 .05 

0.40 

7, 

1594.81 

1450.41 

5.99 

1 399 . 57 

0.32 

83 

1131.77 

1141.54 

0.86 

1156.57 

0.42 

34 

1411.03 

1441 .61 

2.15 

1419.50 

0.56 

83 

1789.12 

1880.52 

5 . o' 

1797.40 

0.4b 

90 

1470.99 

1493 . 53 

l.li 

1484.00 

0.4' 

95 

1810.59 

1 So-'t . 52 

2.89 

1325.60 

0.72 

9? 

2010.0 

2099.19 

4.43 

2023.88 

0.69 

9q 

2223.55 

2560 . 37 

6 . 33 

2259.60 

0 . 63 

iOo 

1724.'! 

1'47.54 

1.32 

1736.16 

0 .60 

IO4 

2054 . 37 

2109.03 

2.60 

2072.37 

0.88 

106 

2254.36 

2342.81 

5.92 

2274.71 

0.90 

IO3 

2471.59 

2609.26 

5.57 

2493.26 

0.8S 

lOn 

2702. G? 

2907.95 

7 . o2 

2724.24 

0.82 
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Specifically,  an  examination  of  Table  S reveals  that  deviations  between 
the  model  without  centrifugal  distortion  and  experiment  increase  dras- 
tically with  increasing  rotational  energy.  For  example,  the  9o  and  lOj 
states  deviate  from  the  e.xperimentally  derived  values  by  only  1.11’j  and 
1 . 32°j  whereas  for  the  Sg  and  IO^q  states,  the  deviation  is  0.521  and  '.021. 

.N'ext  we  consider  the  predicted  level  structure  of  the  correlated 
bender.  Ive  first  note  that  all  values  tlirough  J = 10  lie  well  within  11 
of  the  experimental  values,  .is  some  of  these  lie  above  the  (OlOJ  J = 0 
vibrational  state,  it  is  clear  that  this  model  presents  an  accurate  pic- 
ture of  the  structure  of  the  molecule.  To  obtain  a more  detailed  account 
of  the  model's  success,  we  again  examine  centrifugal  distortion.  For  the 
Sg  state,  the  correlated  bender  accounts  for  S5.12  cm'*  or  90.91  of  the 
effects  of  centrifugal  distortion.  Similarly,  for  the  10^ 3 state,  it 
accounts  for  ISS.'l  cm' ^ or  S9. 21  of  the  energy.  Tl-.us,  the  correlated 
bender  accounts  for  about  901  of  the  effects  of  centrifugal  distortion 
in  HgO.  In  Table  9 we  summarice  these  considerations  for  the  imnortant 
high  energy  states. 

.An  examination  of  Table  9 reveals  that  the  correlated  bender  accounts 
for  S5  to  90l  of  the  energy  due  to  rotational  distortion  except  for  the 
anomalous  states:  S4,  95,  IO4  and  IO5,  i.e.,  the  relatively  low  t states. 

To  account  for  this  discrepancy,  we  note  t;\at  we  have  not  included  Corio- 
lis forces,  which  lead  to  energy  changes  tl'.at  vary  with  the  rotational 
state.  Furthermore,  this  rotational-vibrational  interaction  enters  into 

t!^e  molecular  Hamiltonian  bv  means  of  terms  of  tie  form  ".4  - . Re- 

a aa  a 

writing  this  is  a somewhat  more  transparent  fashion. 
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wliere  can  be  app I'oxiinatcJ  as  a Jiagonal,  I'isiJ  tensor  for  the  sake  of 


simplieit}’.  Here  anJ  r\  refer  to  the  vibrational  nioJes.  Since  the 
benJini;  mode  involves  the  greatest  deviation,  we  can  take  to  be  paral- 
lei  to  tiic  hydrogen-ludrogen  separation.  Taking  1’^  as  one  of  the  OH 
stretch  modes,  it  immediately  follows  that  t!ie  component  of  J that  is 
parallel  to  the  0-rotation  a.vis  is  of  greatest  importance.  However,  J 
is  parallel  to  tliis  a.xis  only  for  the  low  t states,  and  tlu'refore  it  is 
titese  states  rluit  we  describe  most  poorly  with  respect  to  the  Coriolis 
interaction . 

Table  9 

Effects  of  Centrifugal  Distortion 


J, 

Centrifugal  Distortion 

Correlated  Bender 

.Vniount 

bo 

30.97 

:s.  '3 

9: . "0 

'5 

30.09 

25 . 23 

S3 . S5 

77 

55 . bO 

51.04 

91  .SO 

:9.9s 

: : . 19 

■4.0: 

91.40 

S3. 1: 

90.94 

9_- 

53 . "'3 

40.': 

75.-9 

9- 

SO.  19 

"5.31 

S4.40 

99 

1 lo.s: 

1:0."" 

90. 0: 

10:. 

54.00 

5o . 00 

o-.O- 

10, 

SS . 45 

oS.  10 

-".00 

lOe 

13".0" 

II0.09 

S4 . "o 

lOio 

:o5 . So 

1S3.'1 

S9. 74 

rinally,  in  Table  10  we  have  compared  theor\  ai'.d  experiment  for 
rotat Lonal -vibrational  levels  that  lie  beyond  .1  = 10.  An  examination  of 
tins  table  reveals  tr.at  the  model  of  a correlated  bender  .iccounts  for 
these  states  within  1 to  1.5o  .iccur.ic)'.  In  \ iew  of  tlie  fact  tliat  we  are 

:i 


1 


Table  10 


Comparison  of 

iligh-Hneri,/  Rotational 

-Vibrational  Levels  in 

1,000 J State 

J 

I 

ObserveJ 

Correlated  Bender 

°o  Dev. 

lio 

:i45.01 

2I0I .'4 

0.8" 

115 

:5::.4o 

2549.74 

1.08 

119 

29 3S . 5o 

29'3.0' 

1.18 

nil 

5160,97 

52 1 b . t)0 

l.'b 

1-0 

243'. 0: 

24b0 . So 

0.95 

1-2 

1 J 

0 

1 J 

2b42.65 

1.14 

1C4 

:si3. oi 

2848. b" 

1.25 

139 

292~.5S 

29b5. 20 

1.22 

153 

512S.25 

5 1 ” 1 . 5b 

1 . 38 

iSs 

35  IS . 20 

3598.26 

1.50 

15? 

35S4 . 00 

5o40. 19 

1.5" 

1^-6 

2S80.94 

2900.54 

O.fa- 

14.2 

3101. S' 

5155.81 

1.04 

14o 

5 2 1>  0 . 5 0 

5309 . 54 

1.32 

15-3 

5085 . 9 2 

5111.93 

0.91 

dealing  with  a four  parameter  model  that  has  been  fit  to  only  five  pieces 
of  data,  the  demonstrated  accuracy  is  satisfying.  .Note  that  many  of  these 
rotational -V ibrational  states  lie  above  3000  cnrS  which  is  about  equal 
to  the  (030)  state  (ref.  SI. 

We  note  that  the  off-diagonal  matri.x  elements  of  the  rotation  con- 
stants govern  the  coupling  between  the  various  vibrational  states  that 
arise  from  molecular  rotation,  i-'urtliermore , these  quantities  cannot  be 
directly  measured  and  yield  valuable  information  on  the  vibrational  wave 
functions  as  they  vary  sensitively  with  the  hydrogen-hydrogen  motion. 

For  e.xample,  for  • 0 .\  1-  one  has 

<0|.A|1>  = C ' dqvo(q}  V -\q/:~)'2' 

where  •i'o(vi)  is  the  ground  [(010)]  state  vibrational  wave  function  and  C 
is  a constant  that  does  not  concern  us  here.  Note,  that  for  <0|B|1>  one 
has 

<oiB|i>  = c'  dqoo(q)  -i-iU)  (-0) 

Finally,  in  Table  11  we  have  compared  theory  and  e.xperiment  for  the 
(010)  state.  The  overall  accuracy  for  the  correlated  bender  is  on  the 
order  of  1 to  1.5’  (ref.  9). 

D.  Tlie  .Role  of  the  .Nuclear  Masses 

In  this  section  we  e.xamine  the  significance  of  the  nuclear  masses  in 
our  correlated  bender  model  of  water.  This  is  best  done  by  investigating 
the  rotational-vibrational  spectra  and  structure  of  the  symmetric  isotopes 


I 
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Table  11 


Comparison  of  Rotational-Vibrational  Energy  Levels  in  C^lb)  State 


Observed 

Correlated 

bender 

Dev . ‘ 

1, 

45 

.39 

4b. 

.41 

1 

15 

-0 

99 

.05 

100 

.26 

1, 

24 

149 

.05 

150 

. / 1 

1, 

15 

219 

.23 

.25 

1 , 

55 

33 

313 

.01 

51b 

.52 

1 , 

12 

■+0 

323 

.45 

555 

.35 

1 

54 

4 

411 

. o3 

41o 

.-0 

1 

26 

■^4 

535 

.01 

559 

• / / 

0, 

39 

5o 

531 

. 85 

559 

.43 

0. 

63 

53 

557 

.03 

664 

.91 

1, 

19 

Ss 

Sll 

. 65 

316 

.35 

0, 

64 

^0 

637 

,9- 

637, 

.43 

0, 

07 

62 

304 

,63 

814. 

,82 

1 . 

26 

64 

959 

.39 

969, 

•7 

. / J 

0. 

93 

60 

1159 

. 65 

1144, 

,59 

0, 

45 

~0 

975 

.07 

939. 

,32 

1 . 

51 

"3 

1129 

.71 

1145, 

,95 

1 , 

2b 

- „ 

1510 

,34 

1321. 

.93 

0, 

35 

7 7 

1515 

, 45 

1520. 

. 58 

0. 

55 

1177 

.16 

119b, 

,25 

1 . 

62 

S4 

1506 

,69 

1524, 

,39 

1, 

17 

S3 

1956 

.4b 

1941, 

.97 

0, 

23 

^3 

1545 

.Ob 

1570, 

,72 

1 , 

DO 

1952 

. 13 

1944, 

,05 

1 , 

15 

9? 

2157 

.99 

2175, 

,90 

0, 

-4 

93 

2599 

.30 

2406, 

.75 

0, 

29 

10; 

1795 

.08 

1325, 

.12 

i . 

79 

IO4 

21-6 

. 56 

220", 

,93 

1 , 

45 

lOe 

2405 

,21 

2450, 

,62 

1 , 

14 

lOs 

2646 

.41 

2bb7, 

, 14 

0, 

78 

D:;0  and  T2O.  Specifically,  we  present  calculations  of  the  A,  B,  and  C 
rotation  matrices  as  well  as  the  rotational -vibrational  energy  levels. 

By  comparing  these  results  with  H2O  a detailed  picture  of  the  role  of 
the  nuclear  masses  in  the  motion  of  tl\e  two  hydrogen  atoms  can  be 
obtained. 

Before  presenting  these  results,  we  first  note  that  our  effective 
potential  V(q}  depends  on  the  nuclear  mass  in  two  ways:  (1}  through  the 
quantity 

n - ^ V 

^ ■ ■ 2 ^ ^aa 
a 

and  (23  the  correlated  bender  coordinate  q.  The  dependence  of  tiie  u's  on 
the  nuclear  masses  is  obvious  and  need  not  be  discussed  in  detail  here, 
although  we  will  present  an  order  of  magnitude  estimate  of  its  mass- 
dependence  below.  On  the  other  hand,  the  dependence  of  q on  the  nuclear 
mass  is  interesting  and  characteristic  of  the  technique  we  are  using  to 
model  water  and  therefore  must  be  discussed,  he  begi;i  by  noting  that  for 
any  model  of  water,  q corresponds  to  a one-dimensional  path  through  the 
potential  energy  surface  of  the  molecule.  This  path  may  be  either  linear 
or  curved.  For  the  rigid  bender,  q is  merely  the  bending  coordinate 
itself  and  can  be  identified  with  the  bending  mode  a.xis  (see  Fig.  1).  On 
the  other  hand,  in  the  correlated  bender  model  of  water,  q represents  a 
curved  path  such  as  that  displayed  schematical i)' in  Figure  2 . The  degree 
in  which  this  coordinate  departs  from  a straight  line  parallel  to  the 
bending  mode  a.xis  depends  on  the  correlation  between  the  bending  mode 


and  the  OH  stretch  modes. 


Figure  2.  q is  regarded  as  a curvilinear  path  tlirough  the 
potential  energy  surface  of  the  molecule  for  the 
case  of  a correlated  bender. 


To  estimate  the  contribution  of  the  first  effect  we  again  use  the 
rigid  bender  model  as  a test.  Specifically,  we  apply  the  nuclear  poten- 
tial energy  V(q3  generated  in  Section  IIB  for  H2O  and  calculate  the  first 
transition  frequenc}'  of  DrO.  he  find  a deviation  of  -11  cra’^  from  the 
correct  transition  frequency.  Thus,  taking  into  account  U,  the  mass 
dependence  of  the  q coordinate  amounts  to  about  0.5‘’6  of  the  transition 
frequency . 

Our  procedure  then,  is  to  redetermine  the  potential  energy  function 
for  both  D2O  and  T2O  so  as  to  fit  the  first  bending  mode  transition  fre- 
quency. The  results  are  displayed  in  Table  12.  (The  results  for  H2O  are 
included  for  convenience.)  We  note  that  the  first  two  terms,  i.e.,  a2 
and  au,  which  dominate  the  structure  of  the  molecule,  vary  by  1%  from 
the  H2O  to  D2O.  The  remaining  terms  display  a greater  variation;  however, 
they  represent  less  than  O.l^s  of  the  energy  associated  with  the  02  and 
au  terms.  Finally,  all  the  odd  parameters  were  cero  and  \'2  traiisition 
was  fit  exactly. 

.N'ext,  we  examine  tlie  various  rotation  matrices.  As  the  rotation- 
vibration  interaction  is  manifested  in  the  off-diagonal  components  of 
the  rotation  matrices,  we  confine  our  attention  to  these  terms.  Specif- 
icallv,  we  examine  the  various  ratios  A.. /A..,  and  C../C..,  where 

' 11  11"^  11  11'  TT  ’1 


= ;<i|A|j>| 

which  are  displayed  in  Table  15. 
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Table  12 


Potential  Energy  Constants  Display 


Parameter 

H20 

D^O 

ToO 

45510 

ifaOlO 

46590 

au 

-50835 

-32580 

-32147 

as 

-25605 

-20508 

-15002 

aa 

1820 

- 26568 

-2Z75t3 

Table 

Rotation 

15 

.Vatrices 

Parameter 

HoO 

D2O 

T2O 

Ac  I'AflO 

0.22 

0.19 

0. 1” 

A02/ Aqo 

0.09 

0. 0t> 

0.05 

A12/A11 

0 . 55 

0.28 

0.25 

Bci/Bqo 

0.12 

0.10 

0.09 

Boa/Boo 

0.02 

0.01 

0.01 

B12/B1 1 

0.  18 

0.15 

0.15 

1^0  i / o'  0 0 

0.01 

0.01 

0.01 

C;2/Cao 

'0 

'0 

'0 

c 1 2 0 U 

0.02 

0.02 

0.02 

r 


■Vn  examination  of  T;  'c  13  reveals  that  there  is  a systematic  trend, 
in  particular,  the  corresponding  off-diagonal  matrix  elements  of  H2O,  D2O, 
and  T2O  all  decrease  with  increases  in  the  nuclear  masses.  Furthermore, 
these  decreases  are  systematic,  as  each  successive  molecule  decreases 
by  S0°(j,  i . e . , 

CAoi/.\oo)j^q  ' 0-8(Aoi  •■'•00)1)20  ' 0)420' 

Finally,  in  Table  14  we  compare  theory  and  e.xperiment  for  10 
rotational-vibrational  energy  levels  in  the  ground  and  (010)  excited 
vibrational  states.  In  general,  agreement  witii  the  e.xperimentally 
derived  values  is  on  the  order  of  C.b^’o,  which  is  satisfactor)'  in  view 
of  the  fact  the  only  e.xperimental  data  used  in  constructing  this  model 
of  D2O  were  tlie  first  transition  frequency  (i.e.,  V2)  and  the  A and  B 
ground  state  rotation  constants  [ref.  Oj . 

Table  14 


Rotational-Vibrational  Energy  Level  Structure  in  D2O  [cm‘^] 


V 

J, 

Observed 

Theory 

% Dev. 

0 

lOo 

903.19 

911.35 

•^0.40 

0 

I02 

100:. S5 

100”. 65 

+0.47 

0 

104 

1114.35 

1120.34 

*0.55 

1 

104 

1 loo . 35 

11””. 64 

+ 0.92 

0 

106 

1241.93 

1249.16 

+0.5” 

1 

lOe 

1 50” . o7 

151”. IS 

+0.30 

0 

1O3 

1532.69 

1590.97 

*0. 60 

1 

lOe 

1464.2 

14”0.2 

+ 0.41 

0 

lOio 

1535.34 

1545.14 

+ 0 . oO 

1 

lOio 

1656 . “9 

1635.11 

-0.10 
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SECTION  III 


A PREDICTR'E  PHENOMENOLOGICAL  MODEL  OF  WATER  VAPOR 
As  we  have  discussed  in  earlier  sections  of  this  report,  centrifugal 
distortion  dominates  the  rotational-vibrational  level  structure  of  light 
asymmetric  molecules.  For  such  molecular  systems,  only  the  lowest 
rotational  states  can  be  described  within  the  confines  of  a rigid  rotor 
model.  Higher  rotational  states  are  strongl)-  coupled  to  the  vibrational 
degrees  of  freedom,  and  nonrigid  effects  become  an  important  component 
of  the  molecule's  internal  structure.  For  e.xample,  in  water  vapor,  the 
rigid  rotor  appro.ximation  breaks  down  by  J = 4,  and  one  must  include 
centrifugal  distortion  within  .some  well  justified  sclieme. 

From  the  discussion  presented  above,  it  is  clear  that  if  one  wishes 
to  characterize  the  high-resolution  infrared  spectral  data  that  are  avail- 
able for  such  molecular  systems,  a scheme  that  readily  incorporates  non- 
rigid  effects  must  be  employed.  One  such  approach  is  the  Watson  rotational 
Hamiltonian,  which  is  a series  e.xpansion  in  t;ie  various  components  of  the 
molecular  angular  momentum.  This  is  essentially  a perturbative  approach 
that  is  beset  with  a number  of  difficulties  tiuit  render  it  inappropriate 
for  current  Air  Force  needs. 

In  this  section,  we  outline  a new  phenomenological  approach  to 
centrifugal  distortion  in  light  asymmetric  molecules,  particular!)-  H^O. 

This  technique  does  not  involve  a series  expansion  in  powers  of  the  angular 
momentum;  instead  it  is  based  on  the  notion  of  fitting  ti'.e  full  rotational 
matrix  (diagonal  and  off-diagonal  components,'  to  the  spectral  data.  Our 
approach,  which  is  nonperturbat i ve  in  r.ature,  has  the  additional  advantage 
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that  the  phenomenological  parameters  obtained  to  characterize  the  spec- 
trum are  related  to  the  potential  energy  surface  of  the  molecule  in  a 
straightforward  fashion.  Furthermore,  as  we  shall  e.xplicitly  demonstrate, 
our  method  is  predictive  in  the  sense  that  it  will  determine  the  position 
of  rotational  levels  that  lie  far  above  those  used  in  the  fitting  s:heme. 
Finally,  as  we  shall  discuss  in  Section  IlIE,  our  technique  offers  the 
possibility  of  constructing  highly  accurate  potential  energ)'  surfaces 
extending  over  regions  of  many  thousands  of  wave  numbers.  This  in  turn 
will  enable  us  to  construct  spectroscopically  accurate  transition  dipole 
moment  matrix  elements  that  are  required  for  a complete  description  of 
the  absorption  spectrum. 

A.  Theoretical  Foundation 

For  simplicity,  we  shall  confine  ourselves  to  symmetrical  nonlinear 
molecules,  in  particular  H^O,  in  which  tiie  reciprocal  moment  of  inertia 
tensor  is  diagonal;  i.e.,  _ . = 5 . Furthermore,  we  shall  suonose 

that  both  Coriolis  forces  as  well  as  vibrational  angular  momentum  are  of 
little  significance  and  can  be  neglected.  In  that  case,  the  Hamiltonian 


reduces  to 

H = 


/'  .1  C ..  2 P,  “ V 

“ a la  X r X 

a .\ 


CO) 


For  situations  in  which  one  cannot  neglect  the  effects  of  vibrational 


angular  momentum  or  Coriolis  forces,  our  method  can  be  suitably  extended. 
Ae  defer  discussion  on  the  point  until  the  end  of  this  section. 

We  can  best  di.splay  our  approach  b>’  first  outlining  the  usual  tech- 
to  solve  equation  (00).  Let  |k-,...>  and  e^. 


nique  that  one  uses 


be 


the  eigenvectors  and  eigenvalues  when  the  molecule  is  in  tlie  state  J = 0, 


but  with  vibrational  quantum  numbers  ''e  emphasize  that  these 

eigenstates  are  unknown  to  us,  as  is  the  potential  V.  We  shall  assume, 
however,  that  the  vibrational  energy  levels  e,  are  known.  If  the 

K 1 , . . . 

molecule  is  rotating,  centrifugal  distortion  will  give  rise  to  a coupli 
between  the  different  vibrational  eigenstates.  As  a consequence,  the 
rotational- vibrat ional  wave  function  can  be  written  as 

J 


J tM 


J 


,ni. 


. > 


ng 


(:i3 


where  |j!0-!j>  are  the  rigid  syrrmietric  rotor  wave  functions  and  the  e.xpansion 

are  obtained  from  a matri.x  diagonalitation. 


J^.n:,... 


coetricier.ts 


^K,ki,.. 


The  energy  associated  with  the  eigenstate  (21)  is 


,ni , . . . ”ni , . . 


aKK  ^ 

k-,.. 

k i' , . . 


^k',ki , 


sJKlr  ‘■:jK'><ki,...I_  |ki,...> 

' a ' ^ ’ aa ' - 


(22) 


where  the  matri.x  elements  <JKj j fully  tabulated.  Equations 

(21)  and  (22)  represent  a complete  quantum  mechanical  solution  to  the 

Hamiltonian  (20),  and  demand  that  the  intranuclear  potential  energy 

surface  V be  known.  In  general,  this  is  not  the  case.  Instead,  observed 

spectral  data  yield  the  eigenenergies  E,  (and  e ).  We  now 

Jr.n^,...  n^,... 

wish  to  examine  Eqiiation  (22)  as  a means  of  cataloging  spectral  lines  as  well 
as  deducing  a potential  energy  surface  V.  To  this  end,  we  note  that  tre 


only  unknown  quantities  that  appear  in  (22)  are  the  matrix  elements  of 

the  inverse  of  the  moment  of  inertia  tensor,  i.e.,  <ki, . . . [y^^i kf , . . . 

These  matrix  elements  along  with  <JK| | JK'>  appear  in  the  Hamiltonian 

matrix  that  is  diagonalized  to  obtain  the  E.  .On  the  other  hand 

° JTjni,... 

all  matrix  elements  of  the  vibrational  Hamiltonian  H^,  are 

SO  that  only  a knowledge  of  the  vibrational  energy  spacings,  etc.,  is 
required.  Of  greatest  importance  are  the  quantities  <ki, . . . | y^^lkj, . . .>. 
This  suggests  an  interesting  and  practical  approach  to  generating  high- 
resolution  line  spectra  for  molecular  systems.  In  particular,  we  propose 
to  use  the  observed  spectral  line  data  to  obtain  the  diagonal  and  off- 
diagonal  matrix  elements  of  u by  a least  squares  fitting  procedure,  .■is 
noted  above,  matrix  diagonal ization  yields  the  required  coefficients 
'^T,ni,  . . . 

. Using  Equation  (22) , we  can  catalog  the  rotational-vibrational 

''K,ki,... 

energy  levels,  and  in  addition,  predict  m;iny  new  high-lying  states  that 
are  not  used  in  the  fitting  scheme.  This  point  will  be  e.xplicitly  demon- 
strated below.  Before  discussing  simplifying  approximations  that  car.  be 
used  in  the  above  scheme,  we  first  note  that  the  diagonal  matrix  elements 
of  the  y are  simple  related  to  the  various  molecular  rotation  constants 
that  are  often  known.  This  can  reduce  the  number  of  unknown  parameters 
required  for  describing  the  spectrum.  We  note  that  the  off-diagonal  matr 
elements  are  not  directly  observable. 

•Next,  we  determine  the  effect  of  the  various  matrix  elements  of  the 

“ua 


on  a given  rotational-vibrational  band.  This  is  essential  if  wc  are 


to  make  any  s imp  1 L t'y  Lnj;  approximat  loiu; . Suppose  we  are  interested  in  the 
rotational  level  structure  of  a particular  vibrational  band  labeled  b>’  tlie 

>^uantum  numbers  k^ Then,  the  effect  of  the  level  labeled  b\’  k{,... 

is  set  by  the  following  parameter 


n < J(J*-1) 

w 


iNhich  is  equal  to  the  ratio  of  tlie  rotational  enerv’y  divided  b>-  the  vibra- 
tional spaciiqqs  between  tlie  two  levels.  It  follows  then  riiat  near-l\'injj 
vibr.itional  levels  have  a stroni;  influence  on  the  structure  of  the  rota- 
tional positions,  i'hus,  a I'ermi  resonance,  sucli  as  exists  between  the 
tlOd)  .ind  (0_’0)  states  ot"  water  will  be  of  extreme  importance.  On  tlie 
otlier  hand,  dist.int  st.ites  will  liave  .1  nei;lii;ible  effect,  and,  if  desired, 
can  be  i ^4110 red. 

R.  Lomnaiuson  tq  Othe^r  .App roadies 

At  tlie  [iresent  time  there  exists  a number  of  otlier  approaches  for 
calculating;  and  cha  rac  te  ri  c in;c  molecul.ir  spectral  line  dat.i.  Ihese  are 
list  ev.1  b e I ow  . 

1.  Us  iiii;  the  Horn-Oppenhei  mer  approximat  ion.  calculate  the  intra- 
nuclear potential  enerj;\  surface,  l-'ollowinj,  this,  one  then 
solves  tlie  Sclired  inqer  expiation  for  .1  rotat  ins-vibrat  inp  mole- 
cule to  ol'tain  tlie  system's  le\el  structure  ,uid  w.i'.'e  functions. 

I.  Construct  a phenomenoloj;ical  potential  cners.c  surface  from 
experimental  !>■  deduced  vibrational  sp.icinj;s  and  rot.ition 
constants,  l•ollow■.nl;  tins,  one  then  solves  the  Sclirb'di  nj;er 
w.ive  equ.ition  for  tlie  rot  .it  1 on.i  1 - v 1 b r.it  1 ona  1 level  structure 
.ind  ei  cenfunct  ions  . 


3.  The  piienomenological  Watson  Hamiltonian 


H^.  = + A - (J--  - bp  J2) 

-V'  - 

(j.^ji  . j£j_m  . 


is  fitted  to  the  observed  rotational  spectrum  of  water  for  a 
specific  vibrational  band.  The  various  coefficients  A,B,C,dj, 
...P  varv  in  a striking  manner  from  vibrational  state  to 

vibrational  state.  As  a consequence,  one  must  re-do  the  fitting 
procedure  for  different  spectral  regions. 

The  first  approach,  altliough  first  principles  in  nature,  is  limited 
to  accuracies  on  the  order  of  0.1  e\  ; i.e.,  several  hundred  wave 
numbers.  Hence,  it  cannot  be  used  for  liigh-resolution  work.  Purther- 
more,  the  computer  time  for  such  projects  is  trul\'  enormous,  and  in  terms 
of  cost  effectii'eness,  is  totall;.’  inappropriate.  Finally,  such  first  prin- 
ciples studies  take  much  too  long  for  Air  Force  needs. 

We  have  used  the  second  approach  successfully  to  model  the  CO,,  CO, 
and  DF  regions  of  the  water  vapor  spectrum  in  Section  IV.  In  general, 
such  a technique  yields  a predictive,  reliable  model  of  water,  which 


vields  levels  that  are  accurate  to  about  0.1  cm' 


Such  a model  vields 


on 


1 


much  useful  information  on  water  vapor  absorption.  However,  for  ultra- 
high  resolution  studies  that  are  currently  required,  they  must  be  replaced 
by  a somewhat  more  phenomenological  model.  As  we  shall  discuss  in 
Section  HID,  the  present  approach  that  we  will  use  for  HF  can  be  regarded 
as  a first  generation  descendant  of  the  second  technique. 

The  third  approach  is  totally  unsuitable  for  Kirtland's  needs.  The 
reasons  for  this  have  been  discussed  in  a previous  technical  report  (ref.  10'. 
C.  Preliminary  Results  on  the  Bending  Mode  Vibrational  Bands 

In  this  section  we  present  some  preliminary  results  that  illustrate 
our  technique.  These  results  can  be  considerably  improved  with  new  iter- 
ative techniques  that  we  will  incorporate  into  our  programs.  Furthermore, 
this  work  did  not  include  enough  of  the  vibrational  states.  Despite  this, 
our  results  are  impressive.  In  particular,  we  are  concerned  with  the 
ground,  (010),  and  (020)  rotational-vibrational  bands  of  water  vapor. 
Accordingly,  we  have  least  squares  fit  the  various  rotational  matri.x  ele- 
ments to  a total  of  ~2  rotational  levels  that  lie  within  these  bands. 

The  specific  states  used  in  the  fitting  scheme  are  e.xhibited  in  Table  15, 
along  with  the  associated  energies  and  deviations.  In  obtaining  this 
fit,  our  Hamiltonian  matri.x  covered  the  space  spanned  by  the  following 
si.x  vibrational  states:  (000),  (010),  (020),  (030),  (100),  and  (110). 
ive  e.mphasice  th.at  rotational  data  from  or.l>’  the  first  three  of  these 
bands  was  used.  Since  the  rotation  matri.x  elements  are  all  symmetric, 
in  the  sense  that 


< k 1 , . . . ; a I i 


<ki'..  . . iw,_. 


1 


i 
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there  are  65  rotation  constants  tliat  enter  into  the  problem  . As  the  ground 
and  first  excited  diagonal  elements  are  already  known,  this  leaves  57 
parameters  to  be  chosen  to  least  squares  fit  72  rotational  levels.  .As  we 
shall  see,  most  of  these  are  not  significant. 

•An  examination  of  Table  15  reveals  that  our  fitting  procedure  was 
successful,  especially  in  view  of  the  severe  limitations  placed  on  the 
states  used  and  the  fact  that  the  Coriolis  force  was  not  included.  We 
feel  that  this  is  adequate  considering  the  fact  that  we  are  fitting  to 
three  different  bands  with  energies  up  to  5405  cm' . The  fit  could  be 
considerably  improved  (especially  with  regard  to  the  (020)  state)  b>' 
including  the  following  vibrational  states  in  our  Hilbert  space:  (040), 
(200),  (002),  (120).  (.Vote,  it  is  not  necessary  to  have  the  level 
values  for  these  states.) 

Next,  we  demonstrate  that  our  technique  is  predictive  in  the  sense 
that  it  yields  the  correct  rotational  levels  for  states  that  lie  well 
beyond  those  used  for  fitting  purposes.  In  Table  16,  we  have  displayed 
the  predicted  versus  observed  rotational  levels  for  states  that  lie  at 
J = 10  and  higher  in  the  ground,  (010),  and  (020)  vibrational  states. 
Typical  accuracies  range  from  C.b!"  to  l°o  of  the  rotational  energy.  .Note 
that  many  of  the  levels  in  the  ground  vibrational  state  have  rotational 
energies  that  are  greater  than  5000  cm" ^ and  therefore  will  couple  strongly 
to  the  (020)  state.  Furthermore,  the  predicted  results  for  the  (020) 
state  are  excellent.  In  particular,  the  12_5  level  has  a rotational  energy- 
of  2242  cm’*,  and  the  predicted  level  position  deviated  b,v  onl>'  0.2°  of  the 
e.xperimental  value.  To  our  knowledge,  there  does  not  exist  a single  theory 

that  accounts  for  all  of  these  states  as  well  as  that  presented  here. 
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Table  15 


Comparison  of  Observed  and  Calculated 
Rotational-Vibrational  bnergy  Levels 


V 

J 

T 

Obs. 

(cn" 

Cal. 

Dev . 

V 

J 

T 

Obs. 

(cm'M 

Cal. 

Ccnt'^) 

Dev, 

0 

1 

-1 

;3.79 

25.79 

0.00 

1 

7 

•5 

2309.89 

2311. "2 

1.35 

0 

1 

0 

37.  14 

57.15 

0.01 

7 

-5 

2392.33 

2392.79 

0.41 

0 

1 

1 

42.37 

42.37 

0.00 

1 

7 

-1 

24o2.S7 

2460.34 

2.05 

1 

I 

-1 

161S.41 

1613.40 

0.01 

1 

7 

5 

2724.30 

2724.41 

0.11 

1 

1 

0 

lo34 . 93 

1634 . S3 

0.01 

1 

7 

5 

2905.43 

2909.72 

4.29 

1 

1 

1 

1640. 4S 

1640.41 

0.07 

2 

7 

-7 

3753 . 60 

3738.54 

0.06 

2 

1 

-1 

31'5.44 

3175.50 

0.  12 

7 

-5 

3379.  34 

3330.25 

0.91 

0 

3 

-5 

13d. 76 

136.72 

0.C4 

2 

7 

-5 

5967.43 

3969. 34 

1.36 

0 

3 

0 

206.50 

206.55 

0.05 

T 

-1 

4052.33 

4052.01 

0.82 

0 

3 

3 

2S5.41 

235.72 

0.31 

7 

7 

3 

4568.64 

4364.37 

3 . 77 

1 

3 

-3 

1731.92 

1731.92 

0.00 

2 

7 

5 

4573.97 

4573.27 

3.70 

1 

3 

0 

1S13.37 

1315.62 

0.25 

0 

9 

_7 

1079.07 

1077.49 

1.60 

1 

3 

3 

1907.71 

1903.35 

0.S4 

0 

9 

-5 

1201.91 

1201.64 

0.27 

z 

3 

-3 

32S9.24 

5239.61 

0.37 

0 

9 

- 5 

1232.91 

12S2.D3 

0.23 

■> 

3 

0 

3337. 6S 

3334 . 37 

3,27 

0 

9 

1 

1477.29 

1477.36 

0.07 

0 

5 

-5 

525 . 54 

524.90 

0.40 

0 

9 

5 

1651.53 

1632.11 

0.73 

0 

5 

0 

303.97 

504 . 15 

0.13 

0 

9 

5 

2009.33 

2008.79 

1.09 

0 

5 

5 

742.07 

745.74 

1 . 67 

1 

9 

-9 

25i 2 . 37 

2511.41 

0.96 

1 

5 

-5 

1920.70 

1921.14 

0.44 

1 

9 

-5 

2313.40 

2822.55 

4.13 

1 

5 

0 

2126.41 

2125.93 

0.51 

1 

9 

“5 

2904.32 

2904.33 

0.44 

2 

5 

-5 

5473. 93 

54  79. 19 

0.21 

1 

9 

5 

5321.1 

3318.9 

1 

7 

5 

0 

3719.50 

3'’1S.51 

1.01 

1 

9 

0 

5752.33 

3732 . 53 

0.05 

0 

7 

-5 

732.41 

732.42 

0.01 

1 

9 

9 

3994.59 

3933.60 

5.79 

0 

7 

-1 

342.35 

342.29 

0.06 

‘i 

9 

-9 

4063. "0 

4067.44 

1.26 

0 

7 

3 

1059 . S3 

IOl'O  . 33 

1.05 

n 

9 

.7 

4265.14 

4265 . 7o 

0.62 

0 

7 

5 

1216. 19 

1213.13 

1.94 

9 

- 3 

4493.31 

4496.23 

2 . 4" 

0 

7 

•* 

1394.31 

139~.0' 

2.26 

•> 

9 

1 

4734.67 

4'S2.28 

2.39 

1 

7 

_ 7 

2130.68 

21SC.94 

0.26 

9 

3 

4992.14 

4995.10 

2.96 

9 

6 

5483.34 

5431.96 

1 .38 
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Table  16 

Comparison  of  High-Lying  Levels  to  Experiment 


V 

J 

Obs. 

Ccm'l) 

Cal . 

Dev. 

(ca'M 

V. 

J 

T 

Obs. 

Cal. 

(c.u* 

Dev . 
(era'-) 

0 

10 

-6 

1437.19 

1437. IS 

0.78 

0 

14 

'y 

3-165.4 

3457  .Ot) 

S.34 

0 

10 

.2 

1616.45 

1615.73 

0.72 

0 

15 

-8 

5083.92 

3077.23 

6.59 

0 

10 

2 

1S7S.45 

1874.94 

0.51 

1 

10 

-6 

3053 . 6 

3065 . 49 

6.89 

0 

10 

6 

2254.36 

2232.25 

2.12 

1 

10 

-2 

3255.91 

3250.65 

3.26 

0 

10 

10 

2702.09 

2632.14 

19.95 

1 

10 

2 

3565.3 

35o0 .52 

4.78 

0 

11 

-5 

1313.47 

1812.64 

0.83 

1 

10 

6 

3997.3 

3993.2’ 

4 . 53 

0 

11 

-3 

1399.21 

1898.01 

1.20 

1 

10 

8 

4211.0 

4230.98 

10.02 

0 

11 

1 

2144.46 

2142.93 

1.48 

1 

11 

-5 

3323. 55 

3322. 38 

8.83 

c 

11 

5 

2522.46 

2519.04 

3.42 

1 

11 

-3 

3437.59 

3491.54 

4.05 

0 

12 

-6 

2105.87 

2104.76 

1.11 

1 

11 

1 

3660.2 

3654 . 9 

5 . 3 

0 

12 

2300.67 

2298.61 

2.06 

1 

11 

5 

4266.05 

4257.96 

3.09 

0 

12 

2 

2612.94 

2610.54 

2.40 

2 

10 

-10 

4260.36 

4256.98 

3 . 3S 

0 

12 

4 

2813.61 

2803.95 

4.66 

2 

10 

-8 

4480.39 

4480.37 

0.43 

0 

13 

-7 

2415.95 

2412.62 

1 . 55 

2 

10 

-4 

4752.74 

4755.71 

2.97 

0 

13 

-J. 

2586.5 

2584.09 

2.41 

2 

10 

0 

5030.04 

5053. 44 

3.40 

0 

1 5 

*0 

2927.38 

2925.  S3- 

5 . 55 

10 

5237.79 

5240.42 

2 . 65 

0 

13 

3 

5128.25 

3121.70 

6.53 

11 

9 

4714.82 

4-*5.96 

1.14 

0 

13 

5 

3348.2 

3336. 15 

12.05 

•> 

11 

.7 

4905.64 

4909.35 

5.'1 

0 

14 

• 1 

2745.5 

2741.67 

3.33 

2 

11 

-5 

5034 . 39 

5033.95 

4 . 36 

0 

14 

-3 

3085. 0 

3030.39 

4.41 

11 

-3 

5144.41 

5141.05 

2 .76 

0 

14 

+0 

3266.36 

5260.62 

5.74 

12 

-10 

4966.64 

4965.09 

0.55 

-i 

12 

-5 

55S8  .O'* 

5393.59 

4.62 
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The  least  squares  fitting  of  the  rotational  constant  matrix  elements 
to  obtain  rotational-vibrational  levels  for  ground,  (010)  and  (020) 
states  matrix  diagonalization  automatically  yields  the  energy  level  struc- 
ture for  other  bands.  It  is  of  interest  then  to  compare  these  values  to 
experiment.  In  Table  17,  we  have  displayed  the  predicted  versus  the 
e.xperimental  values  for  the  low-lying  rotational  levels  in  the  (110)  state. 

We  emphasize  that  outside  of  the  vibrational  spacings  of  these  states, 
no  experimental  information  was  used.  Finally,  in  Tables  18  through  20  we  have 
compared  predicted  versus  experimental  frequencies  in  the  CO2,  CO,  and 
DF  regions.  In  Table  21  results  for  the  001  band  are  exhibited.  1 

i 

D . Present  Data  Base,  Vibrational  Bands  of  Significance  | 

for  HF  Transitions,  and  an  Iterative  Data  Link  | 

j 

.■\t  the  present  time  tliere  exists  a nujnber  of  sources  of  data  on  i 

water  vapor  transitions.  By  far,  the  most  accurate  are  those  of  Rao  ■ 

J 

(ref.  and  Gordy’s  group  (ref.  b) . Other  data  bases  include  Hall  and  ’ 

Dowling  (ref.  9),  and  Benedict  (ref.  1 1';  . Rao's  results  include  rota- 
tional states  through  J = 15  for  the  following  vibrational  bands:  (000), 

(920),  (100),  (001),  (011),  (110),  and  (050).  Level  positions  are 

reported  to  an  accuracy  of  10'^  cm" ■ and  are  extensive.  .No  data  were  ' 

reported  with  respect  to  the  (010)  vibrational  band  by  this  group. 

Nevertheless,  level  positions  of  10"-  cm**  are  available. 

HF  transitions  occur  from  2 to  5 _m,  i.e.,  they  span  the  spectrum  from 
about  5000  to  500  cm" ■ . The  following  vibrational  hands  arc  of  interest: 

(1)  ground:  (2)  V2 ; (5)  2v;;  (4)  5v2;  (5)  1 , (6)  ji;  (7)  ♦ vs:  and 

(8)  V;  + >2-  -An  examination  of  Benedict's  ESS.A  tables  reveals  that  rota- 
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Table  17 

Comparison  of  Low-Lying  Rotational-Vibrational 
Levels  in  the  (110)  Bands 


i -L. 

j 

T 

Obs. 

(cra"^ ) 

Pred. 

(cm) 

Dev. 

(cm'l) 

% Dev. 
of  rot 
energy 

1 

110 

1 

-1 

5258.40 

5258.42 

0.02 

0.09 

110 

1 

0 

5274.16 

5274.12 

0.04 

0.1 

no 

1 

1 

5279.67 

5279.58 

0.29 

0.6 

no 

_ T 

5504.01 

5504.47 

0.46 

0.67 

no 

1 

0 

5552.01 

5551 . 70 

0.51 

0.52 

no 

-> 

5579.95 

5580.41 

0.46 

0.52 

Comparison  of  Model  to  E.n 


r 


LINE 

Upper  State 

Lower  State 

(0,1,0) ; 

^-7 

(0,0,0) ; 

«— 4 

CO 

(0,1,0)  ; 

"-2 

(0,0,0) ; 

'"4 

(0,1,0)  ; 

Q 

(0,0,); 

1°3 

(0,1,0) ; 

"o 

(0,0,0) ; 

9, 

(0,1,0) ; 

0 

1 

O 

(0,0,0) ; 

^■°o 

(0,1,0) ; 

^-9 

(0,0,0) ; 

'1 

(0,1,0) ; 

(0,0,0) ; 

^-4 

(0,1,0) ; 

"-8 

(0,0,0) ; 

(0,1,0) ; 

CO 

! 

(0,0,0) ; 

(0,1,0)  ; 

0 

1 

(0,0,0) ; 

11_5 

(0,2,0)  ; 

’-7 

(0,1,0) ; 

«-3 

(0,2,0) ; 

‘'-3 

(0,1,0) ; 

■’3 

(0,2,0)  ; 

"-4 

(0,1,0) ; 

-1 

'0 

(0,0,0) ; 

'“4 

(0,1,0) ; 

'-6 

IS 


for  CO;  Laser  Wavelengths 


Observed 
(.c.-n*  ^ ) 

Calculated 

925.49 

925.68 

928.94 

928.11 

944.33 

944.45 

961.03 

964.5 

980.4 

930.92 

1035.14 

1035.86 

1033.07 

1042.37 

1072.19 

1072. <-l 

1076.68 

1077.57 

1091.24 

1093.67 

1108.33 

1106.73 

1129.87 

1130.31 

1143.54 

1145.86 

1118.56 

1122.39 

Deviation 

0.19 

0.07 

0.07 

2.42 

0.52 

:.C3 

4.30 

0.22 

0.69 

2.43 

0.35 

0.44 

1.82 

3.83 

i 

I 

I 

I 

I 


43 


Table  19 

Comparison  of  Model  to  Experiment  at  CO  Laser  Wavelengths 


LINE 


Upper 

State 

Lower  State 

Observed 

Calculated 

Deviation 

(0,0,1) 

'■  5-6 

(0,1,0) ; 

15-7 

1913.24 

1913.75 

0.51 

(0,2,0) 

' 'o 

(0,1,0) ; 

5o 

1903.87 

1904.18 

0.31 

(0,0,1) 

'•  5-4 

(0,1,0) ; 

5-5 

1906.65 

1909.37 

2.72 

(0,2,0) 

(0,1,0) ; 

S 

1920.05 

1921.38 

0.53 

(0,1,0) 

'•  ^°-2 

(0,0,0) ; 

^^-11 

1926.63 

1929.62 

2.99 

(0,0,1) 

'•  1^-11 

(0,1,0) ; 

^5-12 

1917.4 

1919.41 

2.01 

(0,1,0) 

'•  5-1 

(0,0,0) ; 

5-5 

1946.35 

1947.94 

1.59 

(0,1,0) 

'•  ^3 

(0,0,0) ; 

5l 

1957.01 

1957.08 

0.07 

(0,0,1) 

'•  5i 

(0,1,0) ; 

5o 

1965.59 

1969.43 

3.86 

(0,0,1) 

'■  ^6 

(0,1,0) ; 

55 

1983.71 

1985.75 

2.04 

(0,0,1) 

'•  5-2 

(0,1,0) ; 

5-3 

2004.17 

2005.75 

1.58 

(0,0,1) 

'•  5-5 

(0,1,0) ; 

5-2 

2070.24 

2071.69 

1.45 

(0,2,0) 

'•  ^6 

(0,1,0) ; 

54 

2103.00 

2106.67 

3.59 

(0,0,1) 

i. 

(0,1,0) ; 

■^3 

2122.44 

2123.14 

0.70 

(0,0,1) 

■'  5-2 

(0,1,0) ; 

^-1 

2206.67 

2206.07 

: .0 

(0,1,0) 

'•  53 

(0,0,0) ; 

5-7 

2241.93 

2242.53 

:.6 

(0,0,1) 

'•  5-3 

(0,1,0) ; 

5-4 

2301.06 

2301.63 

0. 57 
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Table  20 


Comparison  of 

Model  to 

Experiment 

in  the  DF  Region 

Upper 

Level 

Lower  Level 

Observed 

Calculated 

Deviation 

(0,0,1) 

(0,1,0) 

2510.50 

2512.83 

2.33 

(0,1,0) 

'•  ^^-10 

(0,0,0) 

; 13 

2555.38 

2550.81 

4.43 

(0,1,0) 

(0,0,0) 

' ^2 

2609.72 

2618.72 

9.0 

(0,0,1) 

; S . 

-‘i 

(0,1,0) 

'•  ^-7 

2544.40 

2547.13 

2.73 

(0,1,0) 

(0,0,0) 

; 9 . 

-0 

2690. 38 

2687.43 

2.9 

(0,1,0) 

; II3 

(0,0,0) 

2744.77 

2746.37 

4.4 

(0,1,0) 

; 1)3 

(0,0,0) 

2804.45 

2812.45 

8.0 

(0,1,0) 

: 9, 

(0,0,0) 

; B . 

2543.75 

25:5.76 

2.1 

(0,0,1) 

(0,1,0) 

'•  ^-3 

2497.60 

24?t3. 38 

1.28 

O 

o 

; 10, 

•J 

(0,0,0) 

; 9 

-o 

2690.38 

2604.53 

4.2 

. 45  ' 

! 

i: 

t 
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Table  21 

Rotational-Vibrational  Level  Structure  for  (001}  State 


Observed 

Calculated 

Dev . 

J 

T 

ener.^y  (cra'^) 

energy  (cm  •) 

(cm- " ) 

1 

1 

5796 .98 

5796.98 

0.0 

0 

5849 . 58 

5849 . 56 

0.02 

3 

3 

4050.51 

4050.41 

0.10 

4 

0 

4050.11 

4065 . 78 

0 . 55 

5 

- 3 

4149. S9 

4150.08 

0.19 

5 

3 

4545 . 57 

4545.49 

0.08 

b 

4550 . 70 

4550 . 98 

0.28 

6 

6 

4759.85 

4762.08 

2.25 

7 

-5 

4448.98 

4448.85 

0.15 

- 

1 

4664 . 15 

4665 . IS 

0.02 

- 

5 

4929.07 

4928.44 

0.65 

S 

- 5 

4625.94 

4625.99 

0.05 

S 

0 

4861.80 

4865 . 55 

l.~5 

9 

-9 

4661 .45 

4660.52 

0.91 

9 

0 

5 195 . 4d 

5194.05 

0.59 

9 

6 

5d94.09 

5695 . 28 

0.81 

10 

-6 

5 1 “1 . 05 

51"! .11 

0.06 

10 

0 

5442.09 

544 1.51 

0 . 58 

10 

9 

6335 . 53 

65o5 . 44 

9.31 

11 

-11 

5062 . 00 

50b0 . 52 

1.48 

11 

-9 

5255.25 

5255.02 

2.21 

IVe  envision  the  calculation  of  the  HF  transition  lines  to  proceed 


as  follows; 

(1)  Initial  values  for  the  rotational  constant  matrices  are 
taken  from  our  old  codes. 

(2)  The  rotational-vibrational  spectral  line  structure  is  calcu- 
lated and  compared  to  existing  data  or  new  Kirtland  data. 

(5)  An  iterative  loop  is  set  up  in  which  the  rotational  matrices 
are  varied  until  the  best  values  for  the  spectral  lines  are 
obtained . 

(4}  A Coriolis  matrix  is  included  to  further  refine  the  computed 
spectral  lines. 

(5)  A new  iterative  procedure  is  set  up  until  the  desired 
convergence  is  reached. 

See  Figure  5 for  a flowchart  of  this  procedure. 


E.  Future  A’.mli  cat  ions 


In  this  section,  we  discuss  the  manner  in  which  our  water  vapor 
model  can  be  extended  to  obtain  the  detailed  characteristics  of  the  HF 
absorption  region;  i.e.,  line  shape,  line  shift  and  line  intensity.  We 
first  note  tiiat  the  line  shape  and  line  shift  essentially  arise  from 
rotational  quenching  by  means  of  long  range  multipole  forces.  Hence, 
these  features  of  the  absorption  spectrum  demand  a detailed  knowledge 
of  the  rotational  wave  functions  as  well  as  the  niultipole  forces.  The 
basic  features  of  these  forces  can  be  obtained  from  classical  electro- 
dynamics and  should  not  give  rise  to  any  difficulties.  Furthermore,  the 
rotational  eigenfunctions  are  gei^erated  from  our  computer  package  and 
therefore  pose  no  problem.  In  the  event  that  quantum  dispersion  forces 
are  significant,  a suitable  phenomenology  can  be  developed. 

Line  intensities  demand  a knowledge  of  the  transition  dipole  moment 
matri.x  elements  between  different  vibrational  states.  Thus,  a knowledge 
of  tlie  various  vibrational  wave  functions,  and  therefore  of  the  potential 
energy  surface,  is  required.  To  this  end,  we  note  that  tiie  various 
rotation  constant  matrix  elements  are  sensitive  functions  of  the  vibra- 
tional potential.  For  example,  the  diagonal  and  off-diagonal  matrix  ele- 
ments of  the  A rotational  matrix  are  given  by 

<k:,k^:Aik,',k.'>  = c<kj ^ — i 

1(1  + /2q;  q;‘  2 ■ 

,1  1 i . ' 

skal |k;  > 

COS““ ' 


k-  '> 


where  C is  a known  constant,  is  the  dimensionless  v'ibrational  constant 
for  the  stretch  mode,  and  9 is  the  bond  angle.  The  fact  that  the  A rota- 
tional constant  varies  so  much  from  state  to  state  implies  that  it  depends 
strongly  on  the  vibrational  wave  function  and  therefore  is  a sensitive 
probe  of  the  potential  energy  surface.  Furthermore,  the  large  site  of  off- 
diagonal  raatri.x  elements  such  as  <020jA|l00>  is  a sensitive  probe  of  mode- 
mode interaction  pieces  of  the  potential  energy  surface.  Thus,  by  combining 
the  various  rotational  constants  generated  by  our  fitting  technique  with 
the  J = 0 vibrational  energy  level  spacings,  we  should  obtain  a very 
complete  picture  of  the  potential  energy  surface  and  therefore  of  the 
vibrational  wave  functions. 

Other  molecular  systems  of  interest  tliat  we  feel  can  be  approaclied 
by  means  of  this  ttclinique  are  IIDO,  O3,  .NSO  and  0114. 
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SECTION  IV 


MULITMODE  MODEL  OF  THE  WATER  MOLECULE 
In  this  section,  we  describe  our  multimode  model  of  the  water 
molecule  and  present  results.  In  Section  IVA  we  derive  the  multimode 
Hamiltonian  of  water.  In  Section  IVB  we  describe  the  computer  program 
that  is  used  to  describe  H;0,  and  in  Section  IVC  we  present  results. 

.A . .^lultimode  Hamiltonian  of  W'ate r 

Neglecting  Coriolis  forces  and  vibrational  angular  .momentum,  the 
molecular  Hamiltonian  is  given  by  Equation  (20).  We  now  derive  the  detailed 
form  of  the  kinetic  energy  T = i Since  true  vibrations  must 

be  restricted  to  the  plane  of  the  molecule,  one  need  consider  only  dis- 
placements A.xi,  iyi , ~\2,  '-yz,  and  Ay;  to  describe  them.  Through 

the  use  of  s>anmetry,  a combination  of  the  above  displacements  might  be 
more  convenient.  In  particular,  we  are  looking  for  combinations  that 
are  symmetric  or  antisvmimetric  with  respect  to  reflection  in  the  y-t 
plane  of  symmetry  of  the  molecule.  The  combinations  must  also  preserve 
a fi.xed  center  of  mass,  which  requires  that 


+ My(A.X2  1-X;)  = 0 


(25) 


Mx-yi  ^ ^K-(.-yz  -yz)  = o 

Each  of  the  three  vibrations  illustrated  in  Figure  4 can  be  described  b\' 
a displacement,  Sj^.  In  Figure  4a,  if  both  M^  masses  moved  downward  a 
distance  S;,  must  move  according  to  Equation  *23)  to  preserve  the  center- 
of-mass  position: 


dyi 


M . M . 

- ^ (Ay:  ^ Ay:)  = ^ 2Si 

X X 
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All  Axj^  = 0 in  Figure  4a,  and  the  coordinate  is  symmetric  with  respect 
to  reflection  in  the  y-z  plane  ("flipping"  the  molecule  over  does  not 
disturb  the  direction  of  the  deflections).  In  Figure  4b,  l,X2  = -S2; 

Ax3  = S2,  and  from  Equation  (23; 

M . 

ixi  = - ^ (-$2  + S23  = 0 

X 

all  = 0,  and  S2  is  also  a symmetric  coordinate.  In  the  third  vibra- 
tion, d.X2^3  will  be  some  fraction,  a,  of  S3;  iy2  ^ 3 will  be  a fraction, 
b , such  that 


2X2  = -aSs;  1x3  = -aS3 

4/2  = bS3;  ly3  = -bS3 

Therefore,  Equations  ! 251  yield 


2M 


Ax: 


aS. 


Ay:  = 0 


Since  angular  momentum  about  the  t axis  must  be  zero,  it  can  be  shown 

for  small  vibrations  that  the  M 's  move  along  the  x-v  bonds,  and 

>■  =.  . . 

a . a , , , . 

a = sin  — ; o = cos  — (_4j 


53  is  antisymmetric  with  respect  to  reflection.  \ow  any  arbitrary  vi- 
bration of  the  molecule  can  be  expressed  as  a linear  combination  of  the 
symmetry  displacements,  Sj^: 


2M 


S-  sin  -V 


2M 


_V; 


s- 


-X:'  = 


-S' 


.a 

Sin  ^ 


= -S- 


C* 

S3  cos  — 
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1 


dX3  = S2  - S3  sin 


iy3 


-Si  - S3  cos  -j 


ana  since 


dxi  dijXi  2M, 


-^1  = 


dt 


dt 


S3  sin  ^ 


and  likewise  for  yi,  yi  etc.,  the  vibrational  kinetic  energy  becomes 
2Tvib  = 


= ZMyCeSi^+Sa^+YSj^) 


(26) 


where 


2M  2.^l 

3 = 1+  ; Y = 1 + sin^  4 • 

Common  practice  has  been  to  express  the  potential  energy  function 
of  the  molecule  in  terms  of  bond  length  stretching  and  bond  angle 
bending.  In  Figure  5 ori  and  ir2  are  the  changes  in  the  equilibrium 
bond  length  1.  da  is  the  change  in  the  angle,  a. 


Figure  5.  Example  of  the  changes 
in  the  equilibrium  bond 
length . 
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In  the  limit  of  small  vibrations 


BSi  cos  ^ - S2  sin  y - YS3 


Cl  "X 

gSi  cos  — - S2  sin  y + yS; 


■2eSi  sina/2  - 25?  cosci/2 


Let  us  now  define  a set  of  three  dimensionless  coordinates,  in  terms 
of  iri,  or2,  and  oct,  which  will  serve  to  describe  the  symmetric  stretch, 
the  bending  mode,  and  the  asymmetric  stretch,  respectively: 

ori  + 5r2  eSj  cosa/2  - S2  sin;t/2 


q2  = 'Sa  = 


-23Si  sinv/2  - 2S2  cosct/2 


f r > - in 


Solve  for  in  terms  of  the  q.  : 


— ( qi  cos  — - — q2  sin  — 


qi  sin-4  4 qi  ‘^os  4 


Recalling  the  kinetic  energy,  Equation  ( 2h'i  , wo  calculate 


qi'^  cos^  — - qiq:  cos 


q4  sin“  4 q:q2  sin  4 4 * t q:" 


a . 1 1 • 2 • ^ \ 
- sin  - * j qc  sin  T 1 

1 1 • ‘ a \ 

COS  — + J q;"  I 


and  from  Equation  f2b) 


T = Mni' 


cos^a/2 
3 


+ sin-  -jj  qi 


1 . , 1 • • 1 / sin“3t/2  ’ z j. 

+ y sina  1 - j qiqo  + j ( 3 + ^2  + — 


where 


= NL.2-  (auqi^  + 2ai2qiq2  ■*■  a22q2^  + aasqs^) 


(29) 


ai  1 


1 / 72  - - 

— Icos-  — + p sin 


^f)  = I 


2Mh 

1 . ^ sin- 
0 


a 1 2 = T 5 ina 


0-1) 


2Mo 


sinct 


1 / . 7 a 3 , a \ 1 f . -'‘'h  , 

a22  = jr  (^sin-  3 - a cos-  yj  = 43  ^ - 


cos-  7r 


Y 

4p 


and 


Y = 1 


•■'k  .2a 

- — sin^  TT 


2M 


“h  7 


1 cos-  - 

0 


N'ext,  calculate  the  conjugate  momenta  = 3T/5q^; 

P. 


where 


3 - 


II 

I 

5Io(aiiqi  ai2q2) 

5q2 

3Io(ai2qi  + a22q23 

II 

(50) 


(51) 


0 


and  solve  for  the  in  terms  of  P^: 


qi  = 


q2  = 


I ^22^’!  ■ ^12^2 

I aiip2  - ai2Pl 


qs 


where 


u)  H aiia22  - ai; 


With  E'lquation  (2S;  in  nind,  tor-.T\ 

1 - 2ai2a22?iP2  + aiz^P:- 


qi‘ 


qj* 


qiq: 


, a.:"P2-  - :a;iai2P:P2  * a-^Pr 


I a,;a22PiP2  - ai2a22P."  ' a,ia-2pa^  + ai2^PlP2 


41 . 


Substituting  the  above  into  Equation  !231,  col lect  the  coefficients  of  the 
momentum  type  terms : 

P;^:  a;;a_^2'  * a2  2a:2“  ' -a:2'a2  2 = a;ia2  2''  - ai2“a2  2 

= a2  2- 


P- • a';a'2“  a--a-,t-  - 2ai-ai-2 


a22aii  “ a2ia^2^ 


= aii^. 


P'P5 


2-  "“a;;a-_a22  " •-a-_'a-2a22  * 2a;;a??ai:'  + raj^' 

= - -a^  2^ 


so  that 


a22P;-  - 2a_P;P2  + ai;P2-  P 


S) 


1 


Simplifying  u, 


1 JL  

S 43  \-MqS 


sin-ot 


= iV  Cyy  - X-  sin^a) 

4c“ 


= 4ir[(l  ^ -X  5in"  4)  (l  + 2x  cos-  §)  ' X 


• ? 

^ sin"-a 


1 / . 0 . ^ ct 

= -r^  (l  + 2x  + 4x^  sin-  - cos-  - - x 
4d-  \ « - 


“ sin'^aj 


(1  - :2x)  = t7 


therefore 


= /-  (^22^1^  - -aizPiPz  ^ a:;P^=) 
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This  is  the  final  form  of  the  kinetic  energy  for  a multimode  calculation 
where  the  aj_j  are  given  in  Equation  (50).  Notice  that  it  either  mode  one  or 
two  is  considered  separately,  so  that  no  interaction  e.xists,  ai2  - 


ai;a22  - a-2-  = aiia22 


Therefore , 


p.;  p,:  p, 


0 van  a22  03  3 


1.  Rotation  Constants 


The  inverse  moments  of  inertia  (roation  constants)  are  calculated 


as  follows: 


I = M.X.-  = :h,c- 

yy  i 1 ^ “ 


-•■hI"  y) 


ZMjjr*-  sin"^  = - cos^) 


(54) 


XX 


) M.y.- 
r i'  1 
1 


‘'•.01 


ZM^jb^  + M^Ca  - b)- 


r^(l  + COS0) 


(35) 


where 


M 


mol 


I + I 
x.\  vv 


:m.  * M_ 

II  0 


(56) 


sin(r/Z)  (the  alternate  description  of  the  bending  mode),  then 

•«*  1 ^ 

1 0 


•V- 


XX  M 


mo' 


I-  C - q a/Z) 


1.5") 


I = M.r^ 

W 


In  both  hquatioH'  (54)  and  ^.'3l  aiiu^e . thevai  ;e  r may  var\',  if  the  S'lnmetri 
or  asymmetric  stretch  is  being  con.sidered.  Defiviing  the  average  r"’  as 
the  average  of  the  squares  of  the  separate  bond  lengths, 

r Z * J..Z 

r2  = ^ 

Z(.."+ .ilr;*..Zr:)  .1-;-  - _r  - 

= ^ (5S) 

where 

r;  = * ,:ri 
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Therefore,  Equation  (.58)  yields 


r-  = [x.2  + (q^-  + Zqiqo  + q32) 


V Cqi^  - 2q.q3  + q32)  = + 2qi  + qi^  + q,^) 


In  any  calculation  involving  either  mode  one  or  mode  three,  the  moment 
of  inertia  may  be  factored  into  bending  and  stretching  integrals; 

<ik  ' j.i>  = <i!R(qi^3)  I i>  * <k!(l  - coss)  U> 

where 

R(q;)  = .V.2(l  + 2q.  + q,2) 

R(q3)  = Mh^^CI  + qa-) 

Thus  a two-mode  calculation  with  either  stretch  requires  two  one- 
dimensional integrals.  three-mode  calculation  would  require  a two- 
dimensional  integral  over  R(qi,q3).  Finally, 


A.  . = <i 

ij 


oO 


hi 


j 


B.  Description  of  Computer  Code 

The  niultivibrat ional  .t.oJo  molecular  physics  program  known  as  MIStRE 
is  currently  capable  of  producing  full  three-mode  vibrational  wave 
functions  for  lUO  and  similar  nonlinear  symmetric  triatomic  molecules. 

The  rotation  constants,  from  which  is  calculated  the  rotational  energy 
spectrum,  are  presently  restricted  to  two-mode  analyses,  mainly  for 
convenience,  hliile  several  of  the  routines  comprising  this  program  are 
general  and  may  be  applied  to  much  more  coraple.x  syste.ms,  others  are 
entirely  problem  dependent.  The  modularisation  of  the  code,  however, 
enables  a versatility  that  makes  treatment  of  more  comple.x  models  a 
relatively  easy  task. 

In  the  discussion  that  follows,  the  three  main  products  of  the  cal- 
culation (vibrational  eigenfunctions,  rotation  constants  and  rotational 
energies)  will  be  broken  down  in  terms  of  the  individual  routines  that 
produce  them. 

1.  Vibrational  have  Functions 

The  main  routine  first  calls  SETUP.M,  which  almost  entirely  defines 
the  current  problem.  It  sets  up  the  molecular  geometry  (masses,  angles, 
and  bond  lengths)  , chooses  a coordinate  system,  and  defines  the  kinetic 
energy  operators  in  terms  of  the  chosen  coordinates.  M-\IN'  then  calls 
iVAVEF.N',  wiiich  will  produce  the  eigenfunctions  for  the  molecule  just 
set  up.  WA'.'EF.N  first  needs  the  integrals  that  will  comprise  the 
Hamiltonian  matri.x,  so  it  calls  l.NTGRT  to  calculate  (1)  the  kinetic 
integrals  (?^  t>TC  terms);  (.2)  the  P type  terms  tiiat  will  participate 
in  kinetic  interactions  between  modes;  (5'  the  potential  matri.x  elements 


til 
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(choice  of  power  series  or  Morse  potentials) ; and  (4)  the  matrix  of 
various  coordinate  moments  (<i|q'^!j>)  that  will  describe  the  interactions 
of  the  potentials  between  modes.  These  integrals  are  calculated  and 
stored  for  each  normal  mode.  INTGRT  will  presently  produce  up  to  10 
basis  functions  (Harnonn.;  or  Morse)  per  mode,  but  is  practically  limited 
only  by  storage  and  some  minor  adjustments.  Harmonic  oscillator  matrix 
elements  are  analytically  derived  using  raising  and  lowering  operators. 

Morse  eigenenergies  are  also  analytically  calculated,  but  P tv-pe  terms 
and  coordinate  moments  are  integrated  using  a Gauss-Hermite  12-point 
numerical  integration.  I.\TGRT  is  virtually  problem  independent. 

Having  at  liand  the  integrals  needed  for  the  Hamiltonian  and  the 
parameters  controlling  the  sites  of  the  potential  terms  and  interactions 
(from  input  or  fitting  procedure) , one  needs  only  to  assemble  the  various  j 

terms  into  a coherent  H;imiltonian.  Tliis  is  the  job  of  0RGA.\’2.  ORGAN'!  j 

takes  the  matrix  of  integrals  produced  by  INTGRT  and  stores  them  j 

, 

wherever  required  in  the  multimode  Hamiltonian.  This  routine  is  also 

problem  independent,  and  theoretically  will  handle  any  number  of  modes. 

We  discuss  here  only  fundamental  limitations. 

All  that  remains  of  the  vibrational  calculation  is  to  diagonalite 

the  Hamiltonian.  We  use  a variable  threshold  Jacobi  method,  JACV.\T,  | 

which  gives  reasonable  fast  and  accurate  eigenvalues  and  vectors  for  I 

1 

medium  site  matrices.  The  vibrational  energy  spectrum  that  has  now 
been  produced  may  be  compared  to  experiment  and,  if  desired,  an  iteration 
may  be  set  up  to  fit  the  calculated  energy  spacings  to  experimental 
values.  If  so  instructed,  WAVEFN  will  call  a nonlinear  function  minimi- 

o2 
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zation  routine,  STEPIi',  which  will  fit  the  energy  spacings  to  measured 
values  by  simultaneously  varying  the  potential  parameters.  STEPIT 
is  a standard  routine  using  a direct  search  method  (no  derivatives)  with 
automatic  step  site  adjustment  and  acceleration.  .\lso  adjusted  during 
tl\e  fit  procedure  is  the  frequencN-  of  the  ground  state  vibration.  By 
rainimiting  tiie  ground  state  Hamiltonian,  one  in  effect  optimices  the 
choice  of  basis  sot. 

2.  Rotation  Constants 

If  desired,  M/\1N  will  next  call  ROTMM,  which  calculates  single  or 
double  mode  rotation  constants  (the  bending  mode  alone,  or  in  conjunc- 
tion with  either  stretch  mode).  The  convenience  of  the  present  setup 
is  duo  to  the  fact  tliat  the  double  mode  calculation  can  be  factored 
into  two  one-dimensional  integrations.  A three-mode  analysis  would 
require  a two-dimensional  integral  over  the  stretcii  coordinates.  KOTMAT 
tirst  calls  l.NTKOT,  a separate  entry  point  in  INTGRT,  to  get  the  in- 
verse inertia  moment  integrals.  The  integrals  defined  in  I.XTCRT 
are  problem  dependent  (on  geometry  and  coordinate  system) , but  other- 
wise the  12-point  ilauss -Hcrmi te  integrations  are  general.  Morse  or 
Harmonic  basis  sets  arc  available.  ROTM\T  then  assembles  the  integrals 
into  rotation  matrices,  reduces  tlie  matri.x  to  contain  bands  of  interest 
to  subsequent  calculations,  and  similarity  transform.s  the  reduced 
matrix  to  the  eigenvector  b.isis  found  in  WAA'liHN . 

At  this  point,  MAIN  will  optionally  ,idjust  tlie  molecular  geometry 
to  fit  selected  rotation  constants.  At  present , STI.PIT  varies  the  bond 
angle  and  equilibrium  0-Ii  bond  length  to  fit  exactly  the  experimental 
ground  state  A and  B rotation  constants. 
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3.  Rotational  Energies 

If  rotation  constants  have  been  calculated,  M\IN  will  optionally 
use  them  to  derive  the  rotational  energy  spectrum  for  the  vibrational 
bands  specified  in  the  previous  step.  ROTENG  will  call  IVAXG,  which 
forms  the  vibrational-rotational  Hamiltonian.  KANG  utilizes  the  K.ANG 
transformation,  wliich  factors  the  Hamiltonian  into  four  smaller  matrices 
that  can  be  diagonalized  separately  by  JACV.\T.  This  saves  much  time 
and  space,  but  restricts  the  analysis  to  asymmetric  rotor  type  molecules. 

C . Molecular  Potential  Surface  and  Geometr>' 

It  was  decided  that  the  following  form  of  the  potential  will  pro- 
duce wavefunctions  that  will  best  reproduce  the  rotational-vibrational 
spectrum  of  HiO; 

V(qi.qz)  = D,[l  - exp(-£;qi)]2  + D2[l  - e.\p (-.12^:)  1 " 

* Pi;qi[l  - V'xp(-32qa)]“  + D;)q..ll  - expl,-iiq;)  ]- 

This  potential  involves  six  free  parameters  that  were  least-squares 
fit  to  four  fundamental  (100 , 200 .010 ,020}  and  two  combination  (110,120) 
vibrational  bands.  In  addition,  the  geometric  parameters  of  the  mole- 
cule rQ  and  vq  were  varied  sel f-cons istent ly  during  the  vibrational  fit, 
so  that  the  ground  state  A and  S constants  matched  the  experiment.  Tlie 
static  geometry  of  the  .iiolecule  is 

= 0.9S()2  .A  = 104.53'’ 

The  vibrational  spacings  were  considered  fit  when  within  1 cm"^  of 
e.xperiment.  The  final  values  of  tiie  potential  parameters  are  (in  cm"“}; 
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Dl 

= 65284 

qi 

= 2.46 

Di2  = 

-24559 

D, 

= 56073 

q2 

= 0.7095 

D21  = 

- 5326 

A stability  analysis  was  performed  to  see  the  effect  of  expanding 
the  original  basis  set  of  three  Morse  functions  per  mode.  It  was 
apparent  that  no  significant  difference  occurred  that  would  warrant 
a larger  set,  and  indeed,  the  predicted  050  vibrational  spacing  was 
within  0.0t)'’j  of  e.xperiinent . 

Table  22  shows  101  rotational  energies,  calculated  from  rotation 
constants  predicted  by  the  above  potential. 
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Rotational  Energies  Calculated  from  Rotation  Constants 
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Table  22. 
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